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1. INTRODUCTION 


THE first step in any attempt at a theoretical evaluation of the specific heats 
of crystals is to determine the nature of the disturbance in their structure 
produced by the thermal agitation. In an earlier publication! and more 
fully again in the address to the Conference of Nobel Laureates held at 
Lindau in June 1956 which was printed in these Proceedings,” this problem 
was discussed. A general theory of the specific heats of crystals embodying 
the ideas expressed in those memoirs has since been formulated and pub- 
lished. The theory leads to a determination of the spectroscopic behaviour 
of crystals and simultaneously with it to an expression for their specific heats 
as a function of the temperature. The success of the theory in achieving 
these aims has been demonstrated by a series of fully worked-out examples 
and a comparison with the facts of experiment.* 


The principles on which the present approach to specific heat theory 
is based may be briefly sketched here. While the classical mechanics is 
made use of in the theory, stress is laid on the necessity of paying due regard 
at the same time to the basic principles of thermodynamics and to the funda- 
mental notions of the quantum theory. With these requirements in view, 
the crystal is regarded as an assembly in thermodynamic equilibrium of sets 
of similar oscillators in great number, the energies of vibration of which are 
individually quantised. The identification and enumeration of these oscilla- 
tors and the determination of their modes and frequencies of vibration is 
accomplished by taking note of the three-dimensional periodicity of struc- 
ture characteristic of crystals and applying to the elements of that structure 
essentially the same methods as those which are so successful in the field of 
molecular spectroscopy. The determination of the thermal energy of the 
crystal as a function of the temperature is then effected by the same pro- 
cedure as that indicated by Einstein in his fundamental paper® of 1907 intro- 
ducing the quantum theory of specific heats. 
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The present approach to specific heat theory rejects completely the 
ideas regarding the subject embodied in Debye’s well-known paper® and 
also those set out by Born and his school in numerous publications.” As 
the present generation of physicists and chemists has been brought up to 
believe in these latter theories, it would not be superfluous to explain in some 
detail why their rejection is necessary. That indeed is the purpose of the 
present memoir, which is sought to be achieved by a searching examination 
of the premises underlying the theories of Debye and Born and of the conse- 
quences to which those theories lead. It emerges from the examination 
that those premises are untenable and that the conclusions derived from 
them are false. 


2. THE HYPOTHESIS OF JEANS 


The parent of the specific heat theories of Debye and Born and also 
of other similar theories was a publication® by J. H. Jeans which appeared 
in the Philosophical Magazine in the year 1909. Though the title of that 
paper indicated temperature radiation as its subject, actually it was devoted 
for the most part to a discussion of the thermal behaviour of material bodies 
on the basis of classical mechanics. It put forward and sought to establish 
a proposition which may be stated as follows :—the thermal energy of material 
bodies is identifiable with the sum of the potential and kinetic energies of regular 
trains of waves traversing them and forming stationary wave-patterns filling 
their volume. The reader familiar with Debye’s specific heat theory will 
immediately recognise this as the thesis which Debye accepted and elabo- 
rated in his attempt to explain the experimental results for the specific heats 
of elementary solids at low temperatures. It is not surprising in these cir- 
cumstances that Jeans became an enthusiastic supporter of Debye’s theory 
and by his advocacy helped to promote a general belief in the validity of the 
thesis on which that theory was based. 


The waves considered by Debye in his paper were the longitudinal and 
transversal waves in an isotropic solid contemplated by the classical theory 
of elasticity. The formule set out in the paper were based on an analysis 
of the modes of vibration of an elastic solid sphere whose external surface 
was assumed to be rigidly fixed, followed by their enumeration in terms of 
their frequencies of vibration. Debye restricted himself for the most part 
to a consideration of the cases of monatomic solids exhibiting cubic symmetry. 
No such restriction is contemplated in the theory put forward and developed 
by Born and his school, which claims to be valid for all crystals. While 
Debye considered it to be adequate for his purpose to regard the vibrating 
solid as a continuum, the Born theory claims to take account of the discrete 
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lattice structure of the crystal. The waves considered in it are also of a more 
general kind than those recognised by Debye. Nevertheless, the theory of 
Born, like that of Debye, is based on an acceptance of the same thesis, viz., 
the hypothesis of Jeans. 


3. THE FALLACY OF THE HYPOTHESIS 


Simple considerations of a general character suffice to show that the 
hypothesis of Jeans is a wholly misconceived idea. We have only to remark 
that the specific heat of a solid is an atomistic property, since its value at any 
given temperature is completely determined by two quantities, one of which 
is the number of atoms comprised in unit mass of the substance, and the 
other is the effective average frequency of vibration of the atoms at that 
temperature. Any attempt to explain this situation theoretically has, of 
necessity, to base itself on an atomistic approach to the problem. In other 
words, we have to consider the individual atoms in the crystal and investi- 
gate their dynamic behaviour. Per contra, wave-motions and the stationary 
wave-patterns resulting from the presence of an external boundary limiting 
their travel are macroscopic concepts, which, it should be emphasised, do 
not require us to take any notice of the atomicity of structure of the material 
traversed by the waves. In these circumstances, it is obvious that the 
attempt to build a theory of specific heats on the assumption that tne thermal 
agitation is identifiable with wave-motions is an illogical and illusory pro- 
ceeding. 


We may usefully enlarge upon the foregoing remarks. The specific 
heat theories based on the hypothesis of Jeans concern themselves with the 
enumeration of the stationary wave-patterns assumed to exist within the 
volume of the solid. But, since the atomicity of structure of the material 
need not at all be considered in the description and enumeration of the 
wave-patterns, it remains completely outside the problem. The identifica- 
tion of the thermal energy of the solids with the sum of the kinetic and 
potential energies of the wave-patterns is therefore incapable of leading us 
to any determination of their specific heats. The difficulty which arises 
here is one of pure logic and is fundamental. In the theories now under 
discussion, the difficulty is sought to be circumvented by introducing an 
additional assumption, viz., that the total number of wave-patterns is equal 
to thrice the number of atoms comprised in the crystal. While this assump- 
tion may appear plausible, it cannot be pretended that such a fixation of the 
number of wave-patterns is a necessary consequence or concomitant of the 
wave hypothesis. Actually, of course, it is introduced in order “‘to make 
the theory work’’, in other words, in order that the results of the theory 
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might bear some resemblance to the experimental facts. In logic, such a 
proceeding would be described as petitio principii, or begging the question. 


4. THERMODYNAMICS AND THE THEORY OF SPECIFIC HEATS 


It can be readily shown that the notions underlying the theories of Debye 
and Born are irreconcilable with the statistical concept of the nature of thermal 
energy which lies at the very base of the science of thermodynamics. A 
crystal is an assembly of an immense number of individual atoms which 
are capable of moving from their positions of equilibrium. Hence, the state 
of the system can only be described in terms of the values of a very large 
number of independent variables. No doubt, the atoms are held together 
by their mutual interactions to form the ordered assemblage which we call a 
crystal and which can be regarded as a unit in the macroscopic sense. But 
the statistical concept of thermal energy forbids us from regarding the entire 
assembly as a single unit in the thermodynamic sense. We cannot, for 
example, identify the thermal energy of the crystal with the energy of sta- 
tionary wave-patterns of any sort. Indeed, the constancy of amplitude and 
coherence of phase relationship over extended regions of space and large 
periods of time which the formation of such wave-patterns presupposes are 
excluded by reason of the very nature of thermal energy. Per contra, it is 
to be expected that the thermal energy of a solid exhibits fluctuations in 
space and in time which preclude any definitive description of it except for 
domains of space and periods of time which are very small in comparison 
with macroscopic standards. 


5. CLASSICAL MECHANICS AND THE THEORY OF SPECIFIC HEATS 


That the theories of Debye and Born are untenable becomes evident 
also when we examine them from the standpoint of classical mechanics. 
For, specific heat theory starts from the well-known theorem that all the 
possible vibrations of a system of connected particles are superpositions of 
a set of normal modes in each of which the particles oscillate with the same 
frequency and in the same or opposite phases. The particles in the present 
problem are the individual atoms in the crystal and hence we are here con- 
cerned with the determination of the normal modes of vibration of the atoms 
in the sense of the theorem just stated. Instead of considering this question, 
the theories of Debye and Born concern themselves with an entirely different 
and indeed irrelevant subject, viz., the normal modes of vibration of a macro- 
scopic solid as determined by the form of its external boundary. Their 
treatment of specific heat theory thus stands revealed as a purely fanciful 
approach devoid of any validity or significance. 
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A straightforward application of the standard methods of classical 
dynamics enables us to determine the characteristic modes and frequencies 
with which the atoms in a three-dimensionally periodic array can oscillate 
about their positions of equilibrium. It emerges from the investigation 
that these characteristic modes can be enumerated and that they are related 
in a precisely definable manner to the structure and symmetry properties 
of the atomic groupings in the crystal. In the theories of Debye and Born, 
on the other hand, which assume that the modes of vibration are determined 
by the external boundary conditions and further, that the total number of 
wave-patterns is equal to thrice the number of atoms included in the crystal, 
we encounter results of an altogether different character. According to 
these theories, the great majority of the stationary wave-patterns which they 
contemplate have “‘wave-lengths” of the same order of magnitude as the 
lattice-spacings of the crystal. But these wave-lengths are all different and 
represent different frequencies of vibration. As a consequence, the “cells” 
of the wave-pattern in which the phase of vibration alternates bear no rela- 
tion whatever to the structure of the crystal, while the wave-patterns them- 
selves appear in immense numbers and are all different from each other. 
Such a description of the dynamical behaviour of the atoms linked with one 
another in the geometrically ordered structure of a crystal is, on the face 
of it, a fantastic misrepresentation. It is, in effect, a reductio ad absurdum 
of the theories which lead up to it. 


6. THE QUANTUM THEORY OF SPECIFIC HEATS 


The ideas which connect the hypothesis of energy quanta with the theory 
of specific heats were very clearly expounded in Einstein’s classic paper® of 
1907. If it be assumed that radiation is emitted or absorbed in energy- 
quanta proportional to the frequency, it follows as a necessary consequence 
that the mechanical energy of the oscillator which thus emits or absorbs 
would itself diminish or increase respectively by finite steps of the same 
amount. Generalising this result, Einstein postulated that the energy of any 
elementary structure capable of mechanical vibration would likewise obey 
the quantum rule. He then coupled this with the concept of thermodynamic 
probability embodied in Boltzmann’s principle and showed that the rela- 
tive probabilities of the oscillator being in different energy states could be 
determined with the aid of the principle. The next step is the evaluation of 
the average energy of an oscillator in an assembly consisting of a great number 
of such oscillators having a common frequency of vibration. This leads at 
once to the evaluation of the thermal energy of the whole assembly. The 
basic ideas of the quantum theory were further clarified and illumined by 
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Einstein in the famous paper® of 1917 in which he showed that the exchange 
of energy between a radiation field and the oscillators located in it can pro- 
perly be described only in terms of statistical concepts. The probability of 
the oscillator being in any given state appears in the paper as the product 
of two factors, one of which is the thermodynamic probability factor given 
by the Boltzmann formula and the other is the intrinsic statistical weight 
of that state. 


The foregoing remarks are made to emphasise that the quantum theory 
of specific heats rests on thermodynamic-statistical concepts which are 
wholly alien to the determinism of macroscopic physics. In other words, 
the theory of energy-quanta possesses a meaning and significance only in 
relation to the behaviour of the ultimate units of which matter is composed; 
per contra it is wholly out of place in any considerations regarding the beha- 
viour of matter in bulk. For example, we can speak of the energy states of 
rotation or vibration of a molecule of benzene, but it would be nonsensical 
to discuss the problem of quantising the vibrations of a tuning fork or the 
rotations of a fly-wheel. 


It will be clear from the foregoing remarks that Debye and Born fell 
grievously into error when they rejected the basic ideas of Einstein’s theory 
and sought to build a theory of the specific heats of crystals on the behaviour 
of macroscopic solids. In doing this, not only did they ignore all considera- 
tions of logic as well as the physical principles bearing on the subject, but 
they also sought to impose on the quantum theory a fantastic and altogether 
meaningless interpretation. 


7. THe SPECTROSCOPIC BEHAVIOUR OF CRYSTALS 


The specific heats of crystals and their spectroscopic behaviour are inti- 
mately related properties since they are both determined by the frequencies 
of vibration of the atoms about their positions of equilibrium. This rela- 
tionship is conveyed more precisely by the statement that the elementary 
oscillators in the crystal which are the carriers of its heat energy are also the 
oscillators capable of absorbing, emitting or scattering the radiations incident 
on it. It follows that any valid specific heat theory contains within itself a 
complete description of the vibration spectrum of a crystal: vice versa, a 
theory of the spectroscopic behaviour of crystals which correctly describes 
the facts of observation in that field necessarily contains within itself all the 
data needed for the evaluation of the specific heat as a function of the tempe- 
rature. 
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From the foregoing remarks, it follows that the theories of Debye and 
Born which are based on erroneous assumptions regarding the carriers of 
the thermal energy in crystals also lead to a false picture of their vibration 
spectra. That is actually the case, for the assumptions made in these theo- 
ries lead to the conclusion that the vibration frequencies of a crystal are 
immensely numerous and are all different from each other and hence that 
the vibration spectrum exhibits a continuous band of frequencies stretching 
over the entire range. Per contra, the atomistic theory of specific heats 
shows that the vibration spectrum of a crystal consists of a finite number of 
discrete monochromatic frequencies accompanied by a residual spectrum 
with a relatively small statistical weight which may be considered as conti- 
nuous and which goes down to very low frequencies. 


Spectroscopy is an exact science, and the vibration spectra of crystals 
are accessible to precise study and investigation by several different tech- 
niques and with a great variety of materials. In every case, the experimental 
results support and confirm the results of the atomistic theory in an un- 
equivocal fashion. They are therefore also an objective demonstration of 
the falsity of the Debye and Born theories. 


8. SUMMARY 


It is shown that the hypothesis, originally proposed by Jeans and 
adopted in the theories of Debye and Born, that the thermal energy of a 
solid is identifiable with the energy of wave-motions in its interior is not 
logically sustainable as a basis for a theory of their specific heats. The 
constancy of amplitude and of phase-relationships in extended volumes and 
periods of time demanded by the hypothesis is irreconcilable with the statis- 
tical-thermodynamical concept of the nature of thermal energy, and the 
idea that the theory of specific heats can be treated as a boundary-value 
problem in the physics of macroscopic solids is therefore misconceived. It 
is shown that such an approach is repugnant both to classical mechanics 
and to the principles of the quantum theory. It must therefore give place 
to a purely atomistic approach which gives us a wholly different picture of 
the nature of the thermal agitation in the solid and of its vibration spectrum. 
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LET (91, P1, £23) and (q2, Po, 22) be the two flows of an incompressible fluid 


of uniform density p and kinematic viscosity v, where as usual gis the velocity 
vector, p the pressure and 2 the force potential from which the external 
forces are assumed to be derived. Following Ram Ballabh! we shall call 
these flows as superposable if it is possible to determine a pressure p, + p. +7 


such that Gh a Ge» Pi + Pe + 7, 2, + 2,) is also a solution of Stokes-Navier 
equations with the necessary modification in the initial and boundary condi- 
tions. In the present note we shall consider the two-dimensional flows due 
to a (i) vortex, (ii) spiral-vortex, (iii) vortex-doublet and (iv) a radial flow 
from the point of view of superposability. We may note that in a recent 
paper Prem Prakash? has studied radial flow but we are reconsidering it here 
in order to bring out some important points about it. 


In the first section we have obtained the general condition for superpos- 
ability of two flows and other necessary expressions working through vectors 
in order to make the treatment independent of the choice of co-ordinate sys- 
tem. In particular, we deduce the condition of superposability of axi- 
symmetrical flows and discuss them from the point of view of superposability. 
As far as the plane flows are concerned we have not obtained any new result 
in this section but the discussion of axi-symmetrical flows appears to be new. 
Besides, we have included this deduction here as it is short and elegant and 
we use the condition of superposability in the text in the form deduced here, 


> 
1. Denoting vorticity by w, the equations of continuity and motion for 
the two flows are: 


e > 
div qa, =0 (1.1) 
> > > > 
H+ grad (407 +2 + @)+ 0.x q%=—veurla, (1.2) 
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> 


Ge > 3 Ps > > -_ > 
We first note that condition of integrability of the equation of motion is 


> 
dw 


- os > 
Ga=@-Viqt vie (1.5) 


as will be seen by taking the curl of both sides of the equation of motion. 


Adding (1.1) and (1.3) we find that the equation of continuity 


. > > . 
div (g, + g2) = 9 (1.6) 
for the superposed motion is satisfied. 


Adding (1.2) and (1.4), we have 
>> > > > 
5 (Ma + 9) + grad [41 + a3)? + Pr Pat + (2 + 2)| 
> > > > > > 
+ (a + 2) x (1 + 42) — [ erad(= + 4 - a2) 
> > > > 
+ Oy XG + oy X | 


> oe 
= — vcurl(w, + we) (1.7) 


Hence the two flows are superposable if 


7? >? > > > > 
grad [+ a. do] + os x i + ox q: = 0| (1.8) 
Or 
> > > > > > 
grad m = — p [grad (q, . 42) + we X Gq + o X qa] (1.9) 


The equation (1.9) determines 7. Performing curl operation on (1.8) we have 


> “= > 
curl (w, X q,) + curl (w, X gq.) = 0 (1.10) 
(1.10) gives the necessary condition for superposability of the two flows. 


From (1.10) we conclude that there exists a scalar function x such that 


+ > > > 
We X Gy + w, X Ge = grad x (1.11) 
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so that from (1.8) we have 
a > 
3“ constant — q; . dz — x (1.12) 


In view of (1.11) and (1.12) it is easy to see that (1.10) is also a sufficient 
condition. 


For the self-superposability of the flow q. p, 2) the necessary and suffici- 
ent condition is 


> os 
curl(w x g)=0 (1.13) 
so that 
> ~ —_ 
w X q = grad X (1.14) 
and 
oe > -_ 
= on 2 (1.15) 


We may note the condition of superposability and the expression for 7 
do not contain ¢ explicitly and hence the difference in the consideration of 
steady and non-steady cases arises only through the condition of integra- 
bility (1.5). 


If one flow, say with suffix 1, be irrotational then the condition of super- 
posability of a rotational flow, with suffix 2, on it is 
+ ~ = ~ 
(qi - V) 2 = (We. VG (1.16) 


Considerable simplification results if we restrict our consideration to 
only the two-dimensional flows. In this particular case the above condi- 
tions reduce to 


io > ~ 
(dz. V) #1 + (q - 7) 2 = 9, (1.17) 
for the superposition of rotational flows, 
os ~ 
(q, . V) #2, =0 (1.18) 
for the superposition of a rotational flow on an irrotational flow, and 
7 = 
(q.V)#=0 (1.19) 


for the self-superposition of a rotational flow. 
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We may notice that the two irrotational flows always satisfy the condi- 
tion of superposability. 


We shall now consider the axi-symmetrical viscous flows. We shall 
take the x-axis along the axis of symmetry and as usual denote the perpendi- 
cular distance from the axis by @. The discussion simplifies considerably 
if we work through Stokes’ stream function % such that 


_ _ 1% _ 1 w& 
2 Sre’ 46 3x satis 
and 
See ee oe 
omen Saat ap oe (1.21) 


In terms of the stream function the condition of integrability (1.5) for 
steady motion reduces to 


ay e "eo = s EM, (1.22) 


while the condition of superposability of axi-symmetrical flows with stream 
functions %, and %, takes the form 


py, o-? E*h, Po, D2 Ey) 
—” a )+5 aa )=0, (1.23) 


where J denotes the Jacobian. We may note that in deducing (1.23) we have 
used the equations of continuity for both the flows. 


From (1.23) we find that the condition of self-superposability is 
p, O* By _ 
i's )=0 (1.24) 


while the condition of superposability of a rotational flow, with suffix 2, on 
an irrotational flow, with suffix 1, is 


J Ye — a® (1.25) 


From (1.23) we conclude that the two flows will be superposable if 
(i) they are irrotational, 
(ii) vorticity of each flow is proportional to @, 


(iii) w,/@ is constant along the stream lines of flow 2, while w,/@ is 
constant along the stream lines of flow 1. 
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From (1.24) we see that an axi-symmetrical flow will be self-superposable 
if its vorticity w is of the form 


w = wf (p) (1.26) 
and then the condition of integrability (1.22) reduces to 
E% = 0 (1.27) 


From (1.25) we note that a rotational flow % is superposable on an 
irrotational flow #,, if in general 


wr, = @ f (ts) (1.28) 


and then the condition of integrability becomes 
s(t» SQ) — vEsa sh (1.29) 


We may regard (1.29) as an equation to determine ¥, when #, is prescribed. 


2. Two-dimensional vortex.—Let a two-dimensional vortex of strength 
K be placed at the origin. Then its complex potential is 


W =iK logZ (2.1) 
so that its velocity potential is 
¢ = — K@ (2.2) 


We investigate the general motion which can be superposed on it. The 
condition of superposability (1.18) here becomes 


. dW _ 
i 0 (2.3) 
so that 
we = F (r), (2.4) 


where F is an arbitrary function of r. 


Substituting this in the condition of integrability (1.5) we get 
Ba —elre ti], (2.5) 
where #, is the stream function of the motion with suffix 2, so that 


#, = — VPPH(r) + G(r), (2.6) 
where 


H()=ri+1 
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and G(r) is an arbitrary function of r. 
Now since 
We = Ais 
We have 
” 1 , ” 1 , 
w= —v0[H + H]+[or+ a’ 
From (2.4) w, is function of r only, this is possible, if 
H’ + : H’ _ 0 
and 
F() =G"+1G' 
i.e., if 
H = Alogr+B. 
Hence from (2.7) 


w, = F() =D+C fj exp($ AZ? + BZ) dZ 
and from (2.6) 
#. = O[Alogr + B] + G(r) 
where G(r) is determined by (2.11) and (2.13). 


the centre at the origin. 


tial is given by Milne-Thompson® 
W = (— m+ iK) log Z 
so that its stream function ¥, is 


yf, = — mé+ Klogr 


irrotational flow. Here the condition of superposability is 


devy —_ 
mrs, +K>7=0 


(2.8) 


(2.9) 


(2.10) 


(2.11) 


(2.12) 


(2.13) 


(2.14) 


From (2.13) we conclude that the vorticity in the motion which can be 
superposed on the flow due to a vortex is constant on concentric circles having 


3. Two-dimensional spiral vortex.—Let a two-dimensional spiral vortex, 
i.e., a source and a vortex, be situated at the origin, then the complex poten- 


(3.1) 


We now investigate the general motion which is superposable on this 


(3.2) 





so 


wh 
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so that 


w, = F (7), (3.3) 
where 


r=re-r®, (3.4) 


Also the condition of integrability in terms of the stream function %, 
becomes 


+ oe We — — (1+ %5)(7F Ff’ ; +1) (3.5) 


The complete solution of (3.5) is 


m2 
t= — (1+ Zs) HO) 0+ GO), (3.6) 
Where G(r) is an arbitrary function of r. 


and 
F’” 
H()=tR +1 (3.7) 
using (3.6) we get 


We = A, = — (1 + m) (H" fs =H’) ( de-*°) 
+(1+g)e* Ko 2a +(a"+! c’)] 3.8) 


From (3.3), we is function of + only and in view of (3.8) it is possible 
only if 


H’ + +H’ =0 (3.9) 
and 
m " 
2e r+(c e =G’) = 0 (3.10) 
but then 
w, = 0 (3.11) 


Therefore we conclude that on the flow due to spiral vortex only an irro- 
tational flow can be superposed. 


We may note that when m = 0, (3.8) reduces to (2.9) because then 
T=PF, 
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4. Two-dimensional vortex doublet—We shall now investigate the 
general flow which can be superposed on the flow due to a vortex doublet 
situated at the origin along the axis of y. The velocity potential here is 


=p (4.1) 


so that the condition of superposability of flow (#3, w.) on it is 


IWe We _ 
r—) + tan 6 ~ foe 0 (4.2) 
The complete solution of (4.2) is 
W, = F (7) (4.3) 
where F is an arbitrary function and 
sin 6 
1 (4.4) 


Also the condition of integrability for the flow with suffix 2 is 


— . dpe 
rsin O55 + Fr cos 8 Ht ve, (4.5) 


The complete solution of (4.5) is 


te = — vH(71) = + G(x), (4.6) 
where G(r) is an arbitrary function of + 
and 
H (r) = F’ /F’. (4.7) 
From (4.6) 
1— ——. 
w= A= — CoP a iG (4.8) 





From (4.3) we find that w, is function of 7 only. This is possible if in (4.8) 


H’ = (4.9) 
and 


G’=0 (4.10) 


but then 


(4.11) 





cr repns nk maa 75 eee 
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On Superposable Flows 


so that only irrotational flow can be superposed on the flow due to a vortex- 
doublet. 


5. Radial flow.—with velocity components* 


ara = + £(8) (5.1) 
(492 = 0 (5.2) 

and vorticity 
w= — 5s". (5.3) 


The condition of integrability of the Stokes Navier equation reduces 
to 


2ff' + v(f% +4f)=0 (5.4) 
We shall now determine the most general irrotational flow that can be 


superposed on it. Let the stream function of the irrotational flow be 4, 
then the condition of superposability is 


fare = (5.5) 
The complete solution of (5.5) is 
ip, = F (a) (5.6) 
Now, since %, is the stream function of irrotational motion, 
Ay = 0 
i.€., 
d 3 " 
(C2) + GS) ] F + Gon =0 (5.7) 
or 
F” (4f’ +f)f' (5.8) 





—eo 4(f’)? + (f")?" 


Prem Prakash, working through cartesian co-ordinates, obtains an equa- 
tion corresponding to (5.8) and his satisfying that equation in a particular 
manner amounts to taking 

F”=0 (5.9) 


and 


(5.10) 
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in (5.7). Hence he has investigated only a particular case of the radial 
motion for which the vorticity is harmonic. Further, in this case from (1.5) 


the condition of integrability will be independent of viscosity and (5.4) 
reduces to 





ff' =0 (5.11) 
which means 
f' =0, i.e., f = constant (5.12) 
so that 
on ~ consent. (5.13) 


Hence the radial flow which he considers is due to either a source or a sink 
under his assumptions. 


We may regard (5.8) as an equation in two independent variables @ and 
w, = — 1/r*? f’ (9) and hence it will be satisfied if each side is taken equal to 
a constant k, say, i.e., if 
ws pr = k (5.14) 
and 


4f+f Of 4 


sre +O = _ 
The complete solution of (5.14) is 
ty = F(w) = pay eet +L (5.16) 
where K and L are arbitrary constants. 
On integrating (5.15) we have 
r= [Phe 4] (5.17) 
re () sin (77 + N) (5.18) 
and 
fa (=) f sin® GC +N)d0+P, (5.19) 
where M, N, P are arbitrary constants, and 
> on = a (5.20) 


k+1 
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Now the superposition of an irrotational motion on the radial flow is 
possible only if f obtained in (5.19) satisfies the condition of integrability 
(5.4). Substituting the values of f’, f”, f'Y in it, we get 


pp A= Bie 





) cos ‘6s + N) cosec® 7 + N) = 0. 


Differentiating the above relation and substituting for f we get 


pr+4 + 2at? — 3a = 0 


t= sin (7? +N) 


where 


and 





a=(2) 8v(A — 1) (A— 2) 


This cannot hold for all ¢. Hence we conclude that it is not possible 
to superpose an irrotational motion on the radial flow in general. 


6. In conclusion, we may note that in the cases of the flows due to 
spiral-vortex, vortex-doublet and the radial flow we prove the non-existence 
of the rotational flows in the case of the first two and of irrotational flow in 
the case of the third which can be superposed on them. Such a conclusion, 
however, is not unexpected and it is in the very nature of the inverse 
approach® which we have adopted here. In the case of the flow due to a 
vortex we could find a family of rotational flows which are superposable on 
it. In these flows the isocurls are concentric circles. 


SUMMARY 


In the first section of the present paper we obtain the condition of super- 
posability working through vectors and in particular give explicitly the con- 
dition of superposability of axi-symmetrical flows. This enables us to make 
some general remarks on the possibility of superposition of two axi-symmetri- 
cal flows. The rest of the paper is devoted to the consideration of the possi- 
bility of superposition of general rotational flow on the flows due to a vortex, 
spiral-vortex and vortex-doublet and of irrotational flow on a radial flow in 
two-dimensions. 


In case of the flow due to a vortex, we find that a family of rotational 
flows for which the isocurls are concentric circles is superposable on it. In 
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the case of the remaining three flows we find that the contemplated types of 
flows do not exist. 
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A NEW SYNTHESIS OF DALBERGIN AND METHYL 
DALBERGIN AND A SYNTHESIS OF 
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IN an earlier paper,! 4-phenyl coumarins occurring in nature were consi- 
dered to be related to benzophenone derivatives. This was supported by 
the ready oxidation of the former into the latter. The oxidation of methyl 
dalbergin (I a) producing 2-hydroxy-4: 5-dimethoxy benzophenone (II a) has 
already been described.1 Under the same conditions, 4-phenyl-7-methoxy 
coumarin (Ib) yields 2-hydroxy-4-methoxy benzophenone (IIb). The 
reverse change of a benzophenone into a 4-phenyl coumarin has been used 
in the past in two cases. Limaye and Ramkrishna? boiled 2: 4-dihydroxy 
benzophenone with acetic anhydride and sodium acetate and obtained 
4-phenyl umbelliferone acetate. Similarly using 2: 6-dihydroxy benzophe- 
none, Seshadri and Varadarajan* prepared 4-phenyl-5-hydroxy coumarin 
and its derivatives. This method has now been applied for the synthesis 
of dalbergin (Ic) and methyl dalbergin (I a) using 2: 5-dihydroxy-4-methoxy 
benzophenone (IIc) and 2-hydroxy-4: 5-dimethoxy benzophenone (II a) 
respectively. In exploratory experiments for working out the best condi- 
tions of the reaction 2-hydroxy-4-methoxy benzophenone (II 5b) was em- 
ployed for conversion into 4-phenyl-7-methoxy coumarin (15). The yields 
by this method at their best are not satisfactory as compared with the applica- 
tion of the Pechmann’s reaction using benzoyl acetic ester.! * 5 


However, the benzophenone method is convenient for the preparation 
of 4-phenyl-5-methoxy-6-hydroxy coumarin, which can be called allo- 
dalbergin. As already mentioned Seshadri and Varadarajan* prepared 
4-phenyl-5-hydroxy and 4-phenyl-5-methoxy coumarins by this method. 
The alternative method of preparing these compounds using §-resorcylic 
ester and benzoyl acetic ester along with aluminium chloride in nitrobenzene 
medium did not work satisfactorily. 4-Phenyl-5-methoxy coumarin under- 
goes nuclear oxidation to yield allo-dalbergin whose methyl ether and acetate 
have also been prepared. 


Ahluwalia, Mehta and Seshadri5 reported that partial methylation 
of nor-dalbergin (4-phenyl-6: 7-dihydroxy coumarin) was successful in 
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ais ( Ar\co 878 = OCH a lH 
| ) (6) R=H; R’=OCH, | 
R (c) R=OH; R’=OCH,  R— 
VY» 


CO.C,Hs 
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providing a convenient synthesis of dalbergin (Ic). It has now been found 
that the method of partial demethylation can also be successfully employed. 
When methyl dalbergin (Ia) is heated at 100° with hydrobromic acid in 
acetic acid, a mixture of dalbergin and nor-dalbergin is formed. But by 
using hydriodic acid at 115-20° for half an hour a better yield of dalbergin is 
obtained. The separation of the mixture could be effected by taking advan- 
tage of the property of nor-dalbergin to form a water-soluble complex with 
borax. 


One of the extraordinary properties of dalbergin and methyl dalbergin 
has been the bright colour reactions they give with magnesium and alcoholic 
hydrochloric acid. In this respect there is resemblance to the flavonoids. 
A study has been made here of the behaviour of a number of 4phenyl 
coumarins in this reaction. As will be seen from the results presented below, 
most of the 4-phenyl coumarins give the colour reaction, but it is most 
prominent with 6: 7-disubstituted compounds and least with 5- and 5:7- 
substituted coumarins. The methyl ethers yield brighter colours as com- 
pared with hydroxy compounds. 


The 3-phenyl coumarin derivatives do not give this colour reaction. 
The following methoxy and hydroxy-3-phenyl coumarins have been examined 
in this connection: 7-methoxy, 7-hydroxy, 7: 8-dimethoxy, 7: 8-dihydroxy, 
6: 7-dimethoxy and 6: 7-dihydroxy. It would appear that there is some- 
thing specific in regard to the substitution in the 4-position. Hence, in order 
to examine if a methyl group in this position will serve the same purpose 
as a phenyl group, a number of 4-methyl coumarins have been tested. Their 
reactions are positive though the colours are less intense as compared with 
the 4-phenyl coumarins. The same effect of position of the methoxy and 
hydroxy substituents is found here too. 


In order to examine the combined effect of 3-phenyl along with the 
4-methyl substituent the following 3-phenyl-4-methyl coumarins have been 
examined and they do not give any colour: 7-hydroxy, 6-methoxy, 6-hydroxy 
and 5-hydroxy. 
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The colour reaction with magnesium and alcoholic hydrochloric acid 
(known as the anthocyanin colour reaction) has been largely used in the field 
of flavonoids; generally a deep red or pink colour is observed. The flava- 
none, plathymenin was recorded by King et al.5% to give rise to pink and 
intense blue colour in this reaction. It has not been used before in the case 
of coumarins. From what has been recorded in the previous paragraphs, 
it would appear that this colour reaction is capable of being used for dis- 
tinguishing between differently substituted coumarins. 


EXPERIMENTAL 


4-Phenyl-7-methoxy coumarin.—2-Hydroxy-4-methoxy benzophenone, 
(10 g.),4»®? fused sodium acetate (5 g.) and acetic anhydride (11 c.c.) were 
heated in an oil-bath (165-70°) for 24 hrs., adding 2-S5c.c. of acetic 
anhydride after every six hours. The mixture was poured into water and the 
product (2 g.) crystallised from alcohol, m.p. 116-17°, undepressed by the 
sample prepared by using the Pechmann’s reaction.! 


Permanganate oxidation of 4-phenyl-7-methoxy coumarin.—4-Pheny]l-7- 
methoxy coumarin in acetone was oxidised with potassium permanganate 
added in small quantities. When the product was worked up following 
the procedure described in a recent paper! it gave benzoic acid, oxalic acid 
and 2-hydroxy-4-methoxy benzophenone. 


4-Phenyl-6-hydroxy-7-methoxy coumarin (dalbergin).—2: 5-Dihydroxy-4- 
methoxy benzophenone’ (10 g.) gave acetyl dalbergin (3 g.) which crystallised 
from alcohol, m.p. 157-58°, undepressed by the acetate obtained from 
natural dalbergin.® On deacetylation it yielded dalbergin, m.p. 209-10°, 
undepressed by the natural sample.® 


4-Phenyl-6: 7-dimethoxy coumarin (methyl dalbergin).—2-Hydroxy-4: 5- 
dimethoxy benzophenone’? (10 g.) gave on cyclisation in a similar way methyl 
dalbergin (2 g.),m.p. and mixed m.p. with the natural sample® 144~-45°. 


Partial demethylation of methyl dalbergin: 


(a) With hydrobromic acid—Méethyl dalbergin (2 g.) in glacial acetic 
acid (20 c.c.) was heated on a boiling water-bath with hydrobromic acid 
(66%; 20c.c.) for 1 hr. The solution was poured into ice and the solid 
product filtered. It was digested with a warm (60-70°) solution of 
borax (5 g. in 30c.c. water). The insoluble portion (0-6 g.) crystallised from 
methanol, m.p. 209-10°; mixed m.p. with natural dalbergin® was un- 
depressed. The borax solution was treated with dilute sulphuric acid and the 
product (0-5 g.) obtained crystallised from alcohol, m.p. 267-68° ; mixed m.p. 
with an authentic sample of nor-dalbergin® was undepressed. 
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(b) With hydriodic acid.—Methy| dalbergin (2 g.) in acetic anhydride 
(20 c.c.) was heated with hydriodic acid (20 c.c.) at 115-20° for half an hour. 
The solution was poured into saturated aqueous sodium bisulphite and the 
solid product filtered. It was separated into dalbergin (1-0 g.) and nor- 
dalbergin (0-6 g.) as in the above case. 


4-Phenyl-5-methoxy-6-hydroxy coumarin (Allo-dalbergin).—It was pre- 
pared by the nuclear oxidation of 4-phenyl-5-methoxy coumarin® (2-7 g.) 
in aqueous sodium hydroxide (3 g. in 60c.c. water) with potassium persul- 
phate (4-8 g. in 100 c.c. water) following the procedure described in a recent 
paper.® It (1-5 g.) crystallised from benzene as colourless stout rectangular 
tablets, m.p. 178-79° (Found: C,71-1;H, 4°8; C,gH,,0, requires C, 71-6; 
H, 4-5%). 


The acetate (acetic anhydride and pyridine method) crystallised from 
methanol as colourless stout prisms, m.p. 114-16°. The methyl ether was 
prepared by means of methyl sulphate and anhydrous potassium carbonate 
in acetone solution. It crystallised from methanol as colourless thick rhom- 
bohedral plates, m.p. 134-36°. 


Colour reactions: 


Procedure.—A solution of about 5-10'mg. of the coumarin was dissolved 
in alcohol (5c.c.) and an equal volume of concentrated hydrochloric acid 
added. The solution was treated with magnesium powder in small lots a 
number of times. 


4-Phenyl coumarins.—(1) 4-Phenyl-7-hydroxy coumarin—pink red chang- 
ing to green and finally fading to dirty yellow on standing. (2) 4-Phenyl- 
7-methoxy coumarin—red quickly changing to green and red on standing. 
(3) 4Phenyl-7: 8-dihydroxy coumarin—pink changing to deep red. 
(4) 4-Phenyl-7: 8-dimethoxy coumarin—red changing to brown and again 
deep red; slowly fades to brown yellow on further addition of magnesium 
and hydrochloric acid. (5) 4-Phenyl-6: 7-dihydroxy coumarin (nor-dalber- 
gin}—pink red changing to green and finally deep red on_ standing. 
(6) 4-Phenyl-6: 7-dimethoxy coumarin (methyl-dalbergin)—same as (5). 
(7) 4Phenyl-6-hydroxy-7-methoxy coumarin (dalbergin)—same as (5). 
(8) 4-Phenyl-6-methoxy-7-hydroxy coumarin (iso-dalbergin)—pink red chang- 
ing to green and finally pink. (9) 4Phenyl-5-methoxy coumarin—red, fades 
slowly. (10) 4-Phenyl-5-hydroxy coumarin—orange yellow; fades slowly. 
(11) 4-Phenyl-5-methoxy-6-hydroxy coumarin (allo-dalbergin)—red; fades 
quickly. (12) 4-Phenyl-5 : 6-dimethoxy coumarin—deep red changing to brown 
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red and finally fades. (13) 4-Phenyl-5: 7-dimethoxy coumarin—red fading 
slowly to yellow. (14) 4-Phenyl-5:7-dihydroxy coumarin—orange, fades 
slowly. 


4-Methyl coumarins.—{1) 4-Methyl-7-hydroxy coumarin—pink with slight 
green tinge changing finally to pink. (2) 4-Methyl-7-methoxy coumarin— 
pink changing to blue and finally to red. (3) 4-Methyl-6: 7-dihydroxy cou- 
marin—red changing to pink and finallyred. (4) 4-Methyl-6: 7-dimethoxy 
coumarin—pink changing to blue and finally deep red. (5) 4-Methyl-5: 7- 
dihydroxy coumarin—no colour. (6) 4-Methyl-5: 7-dimethoxy coumarin—no 
colour. (7) 4-Methyl-5: 6-dihydroxy coumarin—no colour. (8) 4-Methyl- 
5: 6-dimethoxy coumarin-red colour. 


SUMMARY 


Though the benzophenone method of synthesis is not quite advantageous 
for dalbergin and methyl dalbergin, it is more convenient for the prepara- 
tion of allo-dalbergin (4-phenyl-5-methoxy-6-hydroxy coumarin) and its 
methyl ether. These compounds and other 4-phenyl coumarins give 
markedly coloured solutions when treated with alcohol, hydrochloric 
acid and magnesium; a number of 4-methyl coumarin derivatives also give 
similar colour reactions. On the other hand 3-phenyl and 3-phenyl-4-methyl 
coumarin derivatives do not respond to this test which may therefore be of 
diagnostic value. 
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STUDIES ON THE DEPENDENCE OF OPTICAL 
ACTIVITY ON CHEMICAL CONSTITUTION 


Part XLVI. The Rotatory Dispersion of d-Camphor-f-sulphonyl- 
nitrophenylamides (o-, m- and p-), d-camphor-f-sulphonyl phenylene diamines 
(o-, m- and p-) and (o-, m- and p-) phenylene-bis-d-camphor-f-sulphonylamides 


By BAWA KARTAR SINGH AND SHIVA MOHAN VERMA 
Received October 30, 1956 


IN continuation of Part XLIII,! we describe the effect of the substitution of 
nitro, amino and d-camphor-f-sulphonylamido (C,gH,;OSO,NH-) radicals 
in place of a hydrogen atom of the phenyl group of d-camphor-f-sulphonyl- 
phenylamide on the rotatory dispersion of the compounds. 


Position Isomerism and Rotatory Power.—We find that Frankland’s 
lever-arm hypothesis* and its electrostatic modification of Rule® do not hold 
for these compounds in any solvent, whereas Cohen’s! rule is obeyed in all 
solvents for camphor-f-sulphonyl-nitrophenylamides (Tables I-III) and 
camphor-f-sulphonylphenylenediamines (Tables [V-VI) but not for phenylene- 
bis-d-camphor-f-sulphonylamides (Tables VII-[X). Whereas the p-, and 
m-isomerides of d-camphor-f-sulphonyl-nitrophenylamide are dextro rotatory, 
the o-isomeride is levo rotatory. 


Influence of Solvent on Rotatory Power.—It is found that the specific 
rotatory power, [a] “{¢, in different solvents, generally increases with the 
decreasing value of the dielectric constant* of the solvent but in the following 
cases a reverse order has been found: d-camphor-f-sulphonyl-o-nitrophenyl- 
amide in chloroform gives [a] *S = —30-5° and in ethyl alcohol [a] % 
= — 36:0°; m-, and p-isomerides of d-camphor-f-sulphonyl-nitrophenyl- 
amide and o- and m-isomerides of phenylene-bis-d-camphor-f-sulphonyl- 
amide show higher rotatory power in methyl alcohol than in ethyl alcohol. 


‘ Characteristic’ Wavelengths derived from Dispersion Constants of 
Drude’s Equation —The compounds now described exhibit ‘ simple’ rotatory 
dispersion in all the solvents as they can be expressed by the one-term Drude’s 
equation. The ‘ characteristic’ wavelengths deduced from the dispersion 
equation vary from 4415A to 1622 A. In solvents of high and medium 
dielectric constants they vary generally from 2394 A to 1622 A (figures in italics 
in Tables I-[X). These high frequency absorption bands correspond to the 








* Solvents in the decreasing order of the dielectric constant: methyl alcohol (31-5); ethy ! 
alcohol (25-8); acetone (21-5); pyridine (12-4); ethyl acetate (6-11); chloroform (5-2); and 
benzene (2-28). 
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general absorption of the saturated molecules. In continuation’® with our 
previous work the authors have noticed that in five cases the rotatory dis- 
persion can be expressed by Biot’s equation and in other cases the low fre- 
quency of the absorption band is replaced by a high frequency in solvents of 
low dielectric constants, namely, ethyl acetate, chloroform and benzene. 
This result in which the low frequency corresponding to the keto group of 
the compound has been replaced by the high frequency corresponding to the 
saturated molecule and brought about by solvents of low dielectric constants 
needs elucidation by direct measurement of ultraviolet absorption of these 
compounds. 


EXPERIMENTAL 


(d- and -d/) Camphor-f-sulphonyl (o-, m- and p-) nitrophenylamides, 
C,)9H,,;OSO,.NHC,H,NO, and (d- and dJ-) camphor-f-sulphonyl (0-, m and p-) 
phenylenediamines, C,)H,,;OSO.NHC,H,NH,. 
Camphor-f-sulphonylnitrophenylamides were prepared by condensing 

camphor-f-sulphonyl chloride (1 mole.) with nitroanilines (2 moles) in presence 

of toluene according to the method of Morgan and Micklethwait.6 Camphor- 

B-sulphonylphenylenediamines were obtained by reducing the nitro compounds 

with iron and acetic acid. The physical properties of these compounds are 

recorded in Table A (1-6). 


(o-, m- and p-) phenylene-bis-camphor-8-sulphonylamides (d- and dl-), 

Cy>H,;0SO,.NH-C,H,NH :SO.C,,.H;;O0 

These compounds were prepared by heating camphor-f-sulphonyl- 
chloride with phenylenediamines in equimolecular proportions on the water- 
bath for two hours. The products thus obtained were washed with water 
and then extracted with alcohol. They were crystallised from dilute aqueous 
alcohol, when colourless crystals were obtained. The physical properties 
of d-, and dl-isomerides of (0-, m- and p-) phenylene-bis-camphor-f-sulphonyl- 
amides are recorded in Table A (7-9). 

The rotatory dispersion measurements were made in a 2-decimeter 
jacketed tube at 35° C. and are given in Tables I-[X. The values of A, in 
these tables are given in A units. 


SUMMARY 


The rotatory dispersion of the compounds (Tables I-IX) have been 
determined in 7 solvents for 13 wavelengths (A = 6708 A to A = 4358 A) 
in the visible region of the spectrum. It is found to be ‘ simple’ as it can be 


expressed by one-term Drude’s equation, [a] *° % = & It is controlled 


Az — Ay?” 
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generally by a low frequency corresponding to the ketonic group of camphor 
and sometimes by a high frequency ‘ characteristic’ of saturated molecules. 
In solvents of low dielectric constants such as ethyl acetate, chloroform and 
benzene the high frequency corresponding to the saturated molecules gene- 
rally replaces the low frequency pertaining to the ketonic group. The effect 
of position isomerism and the influence of the solvent on the rotatory power 
have been also discussed. 
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ABSTRACT 


Recent experiments on high energy electron-photon cascades suggest 
the possibility of the following multiple processes :— 
(i) Electron scattering with pair creation, 
(ii) Electron scattering with the emission of two photons, and 
(iii) Photon scattering with pair creation. 


The cross-section for the first process predicted by Bhabha has been 
recalculated recently by the standard Feynman technique. Using the same 
technique, we here present the calculations for the cross-sections of pro- 
cesses (ii) and (iii). 

INTRODUCTION 


THE analysis of high energy electron-photon showers in cosmic rays suggests 
the necessity of considering third order processes of the following type: 


(i) Electron scattering with pair creation, i.e., ‘ trident’ production, 


(ii) Electron scattering with emission of two photons. The usual 
‘bremsstrahlung’ of Bethe and Heitler (1934) deals with the emission of 
one photon, 


(iii) Photon scattering with pair creation. In the usual treatment, pair 
creation is due to the annihilation of the photon. Here we postulate the 
emergence of a photon after collision. 


Pair production by electrons was predicted by Bhabha as early as 1935 
while the possibility of process (ii) is mentioned in explicit terms by Heitler 
(1954). Recently Murota et al. (1956) have recalculated the cross-section 
for the production of ‘ tridents ’ (pair + scattered electron) by the Feynman- 
Dyson technique. The experimental data of Koshiba and Kaplon (1955) 
seem to confirm the production of electron pairs by electrons. Recently 
Seeman and Glasser (1956) have observed an event in which a helium nucleus 


311 


312 ALLADI RAMAKRISHNAN AND OTHERS 


having an energy of the order of 10° Gev. gives rise to a shower containing 
very many pairs of electrons. The results of Fay (1956)* at Gottingen seem 
to show that the number of pairs is much greater than can be expected from 
the standard cascade theory. The processes (i) to (iii) we have mentioned 
above may explain some of the new data relating to high energy electron 
photon cascades. 


The object of the present paper is to obtain the differential cross-sections 
for processes (ii) and (iii) using the Feynman-Dyson techniques. The stand- 
ard notation of Feynman will be used and explanation will be given only 
when necessary. 


ELECTRON SCATTERING WITH CREATION OF Two PHOTONS 


(a) Matrix elements and their reduction —Throughout this paper, we use 
the natural units 4 = 1 and c = 1 and the following notation: 


a.b =a bt —Azbz—ayby—azbz: The four dimensional scalar product of 
the vectors*** a and b whose compo- 
nents are (at, ay, Ay, az) and (bt, by, 
by, bz) respectively. 

p=E, D): The energy-momentum four vector of 
electron. Subscripts will be used to 


distinguish between different  elec- 
trons. 


g =(o, 9): The energy-momentum four vector of a 
photon. Subscripts will be used to 
distinguish between different pho- 
tons. 


> 
e =(0,e): The spatial vector representing the 
direction of polarization of a photon. 


m: The rest mass of the electron. 


In addition to those we shall use the ‘dagger’ notation of Feynman 
(1953). The daggered operator a** of any four vector a is defined as 


a= atyt — Gx Vx — AyVy — AzYz 
= ApYo — Ayy1 — As¥2 — Ass (1) 


* Private communication. 
** The daggered operator of a four vector a is represented as ¢ by Feynman. For convenience 
of printing we use the corresponding bold face letter. 
*** For convenience ordinary italics are used for four-vectors. No confusion will arise with 
ordinary numbers since in the paper four-vectors occur only in scalar products. 
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where y’s are the well-known Dirac matrices (see for example, Bethe and 
Schweber, 1955). The y’s are characterised by the following relations: 


YuYv + YeY¥p = 0 pA 
=+2 p=v=0 
=-—2 p = vw = I, 2, 3. (2) 


While the matrices yz, yy and yz are anti-Hermitian, yz is Hermitian. 


In this process, a free electron of momentum p, is incident and we obtain 
after its collision with the nucleus, a free electron of momentum p, and two 
photons of momenta q;, qz. The initial and final state electron wave func- 
tions are given by 

fy, = ue Ps" 


ft. = Up eP.7 (3) 


where u,, U, are the four component spinor parts of the free particle wave 
function corresponding to the momenta p, and pz respectively. 


Pity = MU, Pola = MU, 
P1-Py = P2-P2 = m? (4) 


ris a four vector with components (t, x, y, z). The potentials corresponding 
to the emitted photons are given by 


Aw = Cw e~4,.7 


Aoy = Coy, e~‘4,.7 (5) 
where 
€1-91 = O = 2-4, 91-91 = 0 = Go-Go (6) 
The Coulomb potential due to the nucleus of charge Z is given by 
4 Ze 
A=0,A;=¢=->. (7) 


r| 


As we shall be working in momentum representation, it is convenient to define 
V (q) as 


V@ = gap S AM eeratr 
= BF 5 (adr = 1(Q 8 qd (8) 
Q.Q 


> 
where Q is the spatial part of the momentum and 5 the Dirac delta function. 
Ad 
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We shall choose a co-ordinate system in which the nucleus does not move 
and take the direction of propagation of the photon with momentum q, as 
the z-axis and the plane containing the directions of propagation of the pho- 


tons as the. y-z plane. Let the angle between the directions of propagation 
of the photons be @. Thus 


qi = & (Yt — 2) (9) 
G2 = 2 (yt — yz Cos 9 — yy sin 4) (10) 
Pi = Ex vt — Pix Yx — Pry Yy — Piz ¥z (11) 
Po = Es vt — Pox Yx — Pay Yy — Doz Vz (12) 


The photon beam corresponding to (9) can be resolved into two types of 
polarization which will be designated as type A and type B. 


(A) e: = yx (B) a= Yy- 


Similarly the beam corresponding to (10) can be resolved into two types 
of polarization: 


(A’) e. = yx (B’) e, = yz sin 8 — yy cos 8. 


The lowest order Feynman diagrams for the present process are given 
below. 
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Fic. 1. 





Electron scattering with emission of two photons. 
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There are six different Feynman diagrams for the process yielding six 
matrix elements, My, Mu, Mii, Mis, M,, Mu. In the figure, we have 
shown only three and the other three can be obtained by interchanging the 
photons. For convenience we pair the matrix elements corresponding to 
the interchange of photons. The matrix elements are given by 


A 1 1 
= (03 a el ee 
M, + Mi 4nie*v (Q its {7t 5 0 ee e, 
I I 
es er ger see pe (13) 


, a l l 
Mii -t M,, = 4rie®v (Q) Ug {es Pp. + q: omg ec; Pi +Q Roe m Yt 


l l 
ee ee ee ee ee (14) 
° 1 l 
= — ie = = : oie SES 
M, + Mu = — 4mie*v (Q) iis es 5g a Yt ps gs mi! 
l 1 
ee err "= 


We shall demonstrate the reduction of R.H.S. of (13). 


Rationalisation of 
the denominators yields 


” ; ~ (Yt (P2— Q + m) e.(pi — qi + m)e, 
M, + Mu = — 47ie*v (Q) tg ioe 3 Q)? — m®] [(p, — qn? — m?] 


4% (P2 — Q + m) ei (Pi — Ge tm ‘ 
[((p2 — Q)* — m?*] [(pi — qa)? — m?]S * 


To simplify the numerator, we write 


(16) 


tisyt (Pe — Q + m) e2 (pi — Gi + me, my 
= tgyt (P2 — Q + m) e,(— ep, + 2p1-e: — quer + me,) uy, 
= Usyt (Pe — Q + m) ez (2p, — ques) th. (17) 
We next use the relations 
P2— Q= Ph — G— & (18) 
pia = — ap, + 2p,-a (19) 


and eliminate p, from (16) with the help of (4). Reducing the other part in 
a similar manner, we obtain after some simplification 








Pr 
(1: 





of 
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M, + Mu = — 420° (Q) gp © {2b gt (aurea — Pre) 


2p,-e “ 
eh -€y — Py-€y) + 2e,.€3 + er — 


Py-€2 =) + qe Pi-@y | Pr-O —) 
1h Pieh . Pi-N Pi-d2 Pi -2 


qi€iesq, G2@2€1 Gi 20 
sa 271-42 * 271-4 j 


Proceeding in an exactly same manner, we reduce the R.H.S. of (14) and 
(15). The total matrix element M is given by 


M = — 4rie*v (Q) tz {Aryt + Aayt 1dr + Asyt 2g, + Axe, 
+ Ases + Ageiquyt€eqs + Areoqeyteign + Agwoeiqies 





where 
_ __ Pa- 2 Pi-&1 P2- Cy Pi +2 [ Pas ie P2-€ 
i. me aoe Sk. pe Sr Eee 
: Dida Padi Pina Poa t LDa-4i + page 
P2-@: 2p P2-&s ] 2 
— @1.5 F : SS _ 
sta dis P2- Ga? —— *} Q.Q + 2p,.Q 
+ [Bus a; €2-(41 — Pr) “ wat (42 — Pi) 
2 
eile 
ile Q.Q — 2p,.Q 
A I P2-&s Po-€q €2°41 


'“eereel. “oe a eee = Q 
x (Sars — Parts _ Parts) _ (Parts Pres) 
Pi-91 Pi-d2 Pa Ns 271-42 P2-4i =<. 91 Po- Ga 


lhe Hiei 4 ees 1-42) 


.Q + 2p,-Q \po-G P2-Q2 P2-Q2 


ee (- a> a4 vias) 
Q.Q — 2p,.Q Pi-%1 Pi-42 Pi -Qe 


7 bg! 2 P2-4i 271-91 P2- Qe 





P2-@ 4 Pi-&1 ) 
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A, =~ 201 (Pala 1 Paves *3-4) — _Pr-Caer_ 
# ~ Q.Q + 2p;-Q \po-de © Pa-9r ~ Pa-GJ = Pr 91 Pa Ga 
ee) ee 022 i 2 2 2) 2 2 
a = Q.Q + 2p;.Q \po-4 + P2-Q2 et (Pi-91 P1-42) “s 
A, =— wits +} : 
S$  (2p4-s) (Q.Q — 2pe-Q) © (2p2-41) (Q.Q + 2p;.Q) 
ee. Meare 
4(p1-q2) (P2- 1) 
A, =-— Ri + ; 
_~ (271-41) (Q.Q — 2pe-Q) (2p2- G2) (Q.Q + 2p,;.Q) 
en eee 
(471-41)(Pa-42) 
a si 1 
8 (P2-41) (QQ + 2p1-Q) 
a 1 
wins (71-4) (Q.Q — 2p2-Q) 
— 1 
_ (P1- 2) (Q.Q — 2p2-Q) 
1 
wa= (2-42) (Q-Q + 2p,-Q) a 


Using the relations (9) to (12) we can still further reduce the matrix element 
corresponding to various combinations of polarizations and express it finally 
in terms of the y-matrices and the products of y-matrices. Let M,, Moe, M,; 
and M, be the matrix elements corresponding to the polarizations AA’, AB’, 


BA’ and BB’ respectively. After considerable calculation and simplification, 
we obtain 


My = — 4ze*v (Q) tte {Ginyt + Giaye + Gisyy + Gisve + Gisyt¥xvy 
+ Gigytyyyz + Ginveyz¥x + GisyxyyYz} % (23) 

where 

Ay, = Ay — ww (Ag + Az) (1 + cos 9) + wwe (Ag — Ag) — Arg + An 

Gy. = Ag + Ag — weA, — aA, 

Ays = WW (A,g + A,) sin 8 + wwe (Ayg — Ay) sin 8 

yy = WH, (Ag + Az) (1 + cos 9) + wwe (Ag — Ag) 
+ aH (Aig — Ay) cos 8 





1 
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diy = — WA, sin 8 

Ayg = WW (Az — Ag) sin 4 

Ay, = «Ag + w,A, Cos 0 

dy =0 (24) 


Ao, = Ay 
dog = Ag — oAyg — wyW9 (Ag + A, — Ag + Ag) sin 8 
Qo3 = (w.A; — As) cos 6 
Ao = (Ag —_ woAs) sin 6 
dos = W4 We (— Ag + Az) (1 + Cos @) + W Ws (Ag + Ag— Aio— An) cos 6 
Aog = WoAs 
og = wyAg + wW9(— Ag + A, + Ag + Ag — Aro — Ax) Sin 9 
dog = ww. Ag (1 + sin 6) — w,w2A, (1 + cos 6) — ww. (Ag + Ay) cos 0 
+ wy (Arg + Ar) (25) 
Ag, = Ay 
aso = As oa woAs — WWe Ag sin 0 — @,We A, sin ] —_ 1 WoAi9 sin 0 
ms @1WoAy; sin 6 
As3 = Ag — wrAg 
ds, = 0 
Q35 = WeAs sin ff) — Ws Ag ( + cos 0) + W10A, (1 + cos 0) 
+ aya (Ag + Ay — Arg — Ann) 


A356 = wiA,4 
37 =—=_-_ WoAs cos 0 — WyWo Ag sin 0 W1Wo A, sin 7] 
Agg = WW A, (1 + cos 6) — WWos A, (1 + Cos 6) + W1Wo (— Ag <—- Ag 
a Ajo) + W Ws An cos 0 (26) 


Ay, = Ay + wwe Ag (1 + cos 8) + wywe A, (1 + cos 4) + wwe (— Ag 
+ Ag + Aig — Ay) cos 6 


Aye = Ags = Gag = Ayg = O 

Aaa = A; sin 6 — woA; sin tf] — WW Ag (1 a cos 4) —_ Wi Wo A, (1+cos 0) 
+ ww (Ag — Ag) CoS 9 + wa (— Aig + Ay) 

Aan =__ w,A, ae woA,; a W Wo (A, = A, — Ag — AY — Ai 
+ A,,) sin 8. (27) 
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SpuR CALCULATION 


The differential cross-section for any particular polarization is given by 


do; = 2m a - (normalisation factor) - 4. 
1 


D4 = |Mi|* - (density of final states). (28) 


spins 1 spins 2 


where we have summed over the final spin states of the electron and averaged 
over the initial spin states. In this section, we shall evaluate 


=  |Mil*. 
spins 1 spins 2 

It is well known that if 

M = i, Mu, (29) 
then 

EE. |M*| = Sp (ps + m) M(p, + m) M" (30) 

where 

M’ = y, M. (31) 


Now MM, is given by 


Mj; = ainyt + Gisvx + Gisyy + Gave + UisvtYx¥y + UevtY yz 
+ Ginyt¥z¥x + GsYxVyvz: {(32) 
It is easily verified by direct multiplication that 


Pp, M; = lg 1 + ligvtye + lisytyy + liaytye + lisyx¥y + lighy¥z 


+ ligveve + lisytyxyvyyz (33) 
P, Mj = my 1 + migytyx + Misytyy + Miaytyz + MisyxYy 
+ MigyyYz + Mizyz¥x + MigytyxcY yz (34) 


where 
liy = Ex + PoxGig + Pays + Pozia 
lig = Exdig + PoxGiy + Poyis — Pozir 
lig = Exdis — Doris + PoyZir + Pozi 
lig = Exdig + Pox dig — PoyWig + Pozin 
lis = Exdis — PoxGis + Payie + Dozdig 
lig = Eedig + Poxdig — PoyGia + PozQis 
lig = Exdiz + Poxdia + PoyAis — Poze 
lig = Esdig + Poxig + PayGiz + Pozdis (35) 


t We omit the factor — 4arie* v (Q) since it has no role in the calculation of the spur. We 
shall restore the factor when we substitute for |M,|* in equation (28). 
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My = Edi, + Pixie + Pryis + PrzGia 
Mig = E,dig + Pix din — Pry4is + PrzQiz 
Mig = E,dig + PixGis + Piy%i — Piz%g 
Mig = Eydig — Pirin + Piy%e + Pizin 


Mis = — E,dis — Pixdis + PryFie — PizMis 
Mig = — Exdig — Pixdig + Pryia + Piz 
Miz = — Eydig + Piz dig — PiyGs — PizGe 
Mig = — Ejdig — Pig — Pry — Przis (36) 


and I is the unit matrix. 
Now (30) can be written as 


= 3 |Mil|* = Sp|(p, + m) Mj; (Pi + ™) M;'| 


spins 1 spins 2 
= Sp |p. Mj Pp, M;'| + m? Sp |M{M;j| 
+ m Sp |p, M,Mj'| + m Sp | Mjp.Mj'|. (37) 


Using the well-known properties of y-matrices, it can be proved that 
Sp | Mip.M;'| = 0 
Sp |P2M4M;'| = 0 


Sp |poMjp,M;'| = 42 Imi € (J) 
8 
Sp |MjMj’| = 4 Y a,j? €' (J) 
e(j)= +1 forj=1,2,3,4. (jf)=+4+1 forj=1,5, 6,7. 
= — 1 forj = 5, 6, 7, 8. = — 1 forj = 2,3,4,8. (38) 
Thus we obtain 


SS [Myl? = 4 [4me®w (Q)|? S [agri € (i) + mPay2e’ ()]. (39) 


spins1 spins 2 {=i 
DENSITY OF STATES 
The density of states for a 4-particle system can be easily derived as 
—_ps"dpsd2sp,dpd@2,E;E,p;*dQ, §(40) 
(E—E,—E,) p;*—E; {((p — Ps — Pa)-Ps} 


§ We use the subscripts 3, 4, 5 and 6 to denote the four particles instead of 1, 2, 3,4 just to 
avoid any confusion with the use of subscripts 1 and 2 for the initial and final states of the electron 
in the problem. 


D = (22)° 
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Adapting (40) to our problem, we make the substitutions 
E,= 00, Ps =P, Pa =, Ds = 2. 
We then obtain 
D = (22)* p.*dprdQyu,*dw,dQ,, w7dQ,,.. (41) 
Hence do; is given by 
doy = 9 Na a EY (2ny-® pe®dpyeo, *den dQ, .00,5d0. 


x 2 [4ne% (QE [gmye (j) + may?’ (j)). (42) 


If we are not interested in any particular polarization, we can sum over 
all polarizations. This is done by asumming do; over i. 


We shall consider a particular case when the two photons are weak. 
In this case, the cross-section reduces to a very simple form. The matrix 
element (13) can be expressed purely in terms of q, and gp as 








- 1 
M My, = — 4rie* 
i+ Mu ie (Q) unt 5 —g — qm 
l 1 
xX “> —q—m e+e i ——— m*? Uy. (43) 


Rationalising the denominator and neglecting q, and q, in the numerator 
and q,-g2 in the denominator, we obtain 


— — dmie®y (Oy #27t (Px — ™) 
M, oo My = Arie vD (Q) 4p, Ah a qe) 
(Pi — — m) | ey ey A~ m) | 
- ' 44 
” [es Pr r Py- Qe “ay (4) 
Reducing still further by the method now familiar, we obtain 
My + Mu = — 4nie%y (Q) figytty BE OEE (45) 


In a similar manner, we reduce the other two pairs of elements and obtain 





M, = — 4mje? . Paty _ Pa.) Parts — Pra) 46 
, wire Sam Pe: Pi-917 \P2-d2 = P1- 2 aa 


On making the spur calculation, the cross-section can be written as 


2 E Q 
doy = E*[ serge |v (Ql? hy? 4 (m? + EAE, — pr-Pd)| 


















Multiple Processes in Electron-Photon Cascades 


-e*dw,dQuw, Poy ae Pi-ey 2 
x (2:7) ?a, n a | 











q Pare, Prey, 
e*durdQw, {P22 _ Pr-€2 \* 
x 2m)? 4 
) | aon \ 8 wie 


so that the last two factors can be interpreted as the probability of emission 
of the weak photons, if there is elastic scattering from momentum p, to pe 
(see Feynman, 1953). 


SCATTERING OF A PHOTON WITH PAIR CREATION 


In this process a photon of momentum gq, and polarization e, is incident 
: and after its collision with the nucleus, we obtain a photon of momentum q, 
and polarization e, besides a positron-electron pair of momentum p, and p_ 
respectively. The lowest order Feynman diagrams for this process are given 
below. As in the previous section, we have given only three diagrams since 
the other three can be obtained by the exchange of the positions of the 
photons. 
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Fic. 2. Scattering of photon with a pair creation. 








d 
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With respect to the directions of the arrows and without regard to the 
directions of increasing time, these diagrams look exactly like those for the 
process discussed in the previous section (cf. Fig. 1). However when the 
direction of time is taken into account, we notice that 


(1) p, is the momentum of the electron travelling backward in time, i.e., 
Pi = — P+. 

(2) The photon gq, is absorbed instead of being emitted. 

(3) Ps = P_. 


Hence we can obtain the matrix elements corresponding to these diagrams 
from those of process (ii) if only we replace p, by — p,, q, by — q, and p, by p_. 
However the density of final states is different in this case since the particles 
in the final states are now a photon, an electron and a positron. Thus do; the 
differential cross-section for a polarization of i-th type is given by 


2n 


dy = F505 Dea, 


E (27)*p,2dp,.d2, p_*dp_dQ_w.*dQw, 
+ 


x 2 [4re®v (Q)]? Pi [ligmize (J) + m?ajj?' (7) (48) 


where the prime over the square bracket under summation sign indicates 
that we have to make the following substitutions 


— p..for p,; — E, for E, 
p_for pe; E_ for E, 
and 
= for W. 


If we are not interested in a particular state of polarization of the photons, 
we average over the polarization of the incoming photon and sum over 
polarization of the outgoing photon. This is done by summing over i and 
dividing by 2. 


We hope that this process will be of considerable interest in high energy 
electron-photon cascades since it yields a photon, in addition to a pair. In a 
later contribution we hope to discuss these cross-sections under some useful 
approximations. We also propose to include these processes in the usual 
soft cascade theory and study the fluctuations in the number of electron pairs. 


Two of us (S. K. S. and R. V.) are grateful to the Madras University for 
the Government of India Scholarship and the other (N. R. R.) expresses his 
gratitude to the Madras University for a generous stipend. 
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SUMMARY 


A method of inhibiting shrinkage in nuclear emulsions has been 
developed. It consists of replacing the silver halide in the emulsion, 
which is removed during the fixing stage of the processing, by canada 
balsam. This method can be employed with advantage in the study of 
neutron spectra with nuclear emulsions, as it will yield improved energy 
resolution. Another field in which this process of shrinkage inhibition 
can be advantageously employed is the study of gamma ray spectra, with 
heavy water loaded emulsions. Here again, improved energy resolution is 
possible. In fact, this technique can be employed with advantage in all 
nuclear emulsion work, as it will (a) result in overall reduction of dis- 
tortion in emulsions, (b) improve conditions for the measurement of tracks, 
and (c) improve the stability and keeping qualities of processed emulsions. 


INTRODUCTION 


Ir has come to be recognised that the most generally useful technique for 
neutron spectroscopy in the energy region 0-5-15-0 MeV is the nuclear 
emulsion technique. In this method, the trajectories of proton recoils are 
employed to determine the energies of the incident neutrons. In effect one 
needs to measure the range and orientation of the recoil protons. The ulti- 
mate precision with which these quantities can be measured, together with 
the range straggling of the protons, determines the resolution of the method 
and the energy region in which it can be usefully employed. It has not been 
possible hitherto, to extend measurements on recoil protons to energies below 
0-5 MeV and thus extend the measurable energy range of the neutron spec- 
trum. The reasons are twofold. One is the relatively large number of 
background grains resulting from the gamma radiation which is invariably 
present during neutron exposures. As a result, there is considerable straggl- 
ing in range for low energy tracks. Furthermore, measurement of tracks 
containing less than a few grains is almost impossible due to the difficulty 
of distinguishing random background grain patterns from true tracks. 
Secondly, there is the difficulty of determining with any degree of accuracy, 
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the dip angles of tracks which are smaller in length than 5-5 microns (or 
0-5 MeV protons). This arises from the finite depth of focus of the micro- 
scopes used. 


The finite depth of focus can contribute errors in the measurement of 
long range protons as well. Firstly, it is the first few microns of any recoil 
proton track, which gives accurately the direction of the recoil proton relative 
to the incident neutron direction. The direction is thus determined with 
dubious accuracy in those cases where the initial portion of the proton track 
is not favourably situated for the measurement of the dip angle. Such in- 
accuracy in the measurement of a number of tracks results in the broadening 
of the neutron energy peaks finally obtained. Secondly, the total length of a 
track is composed of many independent segments, and the inaccuracies in 
the measurement of individual segments arising again from the finite value 
of the depth of focus of the microscope, all add up to contribute a dispersion 
in the final values of the ranges. 


It will thus be seen that the precision with which the ranges and orienta- 
tions of the recoil protons could be measured, is limited severely by the phe- 
nomenon of shrinkage of the emulsion during processing. This shrinkage 
results from the removal, during the fixing stage of the development, of un- 
sensitized silver halides. The emulsion layer, thus devoid of the silver halide, 
contracts in thickness during the drying process. Such shrinkage leads to 
several types of distortions: 


1. actual contraction perpendicular to the plane of the emulsion, 


2. relative displacement of emulsion layers parallel to the plane of the 
emulsion, and 


3. in cases where foreign substances like wires and capillaries are pre- 
sent in the emulsion, volume distortion by contraction of the emulsion round 
the wire or capillary. Ifa suitable substance could be found to fill the voids 
in the emulsion, so as to bring it back to almost its original thickness, this 
would minimise the distortions in the emulsion, and enable dip angle 
measurements to be made with an accuracy comparable to that with which 
the horizontal angles are measured. 


Another field in which the desirability of shrinkage inhibition has been 
recognised, is the study of gamma ray spectra with heavy water loaded emul- 
sions. Gibson et al.,\»* have shown that this method of using heavy water 
loaded emulsions can be valuable for the measurement of gamma ray energies 
in the interval 4-12 Mev. In this method, the range and orientation of pro- 
tons resulting from the photodisintegration of the deutrons from the heavy 
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water are measured, and are used to infer the energies of the gamma rays. 
Here the precision with which the dip angles of the protons can be measured 
is limited by the shrinkage of the emulsions, to an even greater extent. 
Whereas the shrinkage factor, i.e., the ratio of the thickness of the emulsion 
during exposure to that after processing is ~ 2-5 in neutron spectroscopy, 
the value of this factor in gamma ray spectroscopy lies between 9-5 and 11-0. 
This is because the nuclear plate, at the time of exposure to gamma radiation, 
is soaked in heavy water and is swollen to about four times its initial thick- 
ness. The accuracy attainable in dip angle measurements is therefore even 
more restricted. Gibson et al.? have indicated that even if the final thickness 
of the processed emulsion could be increased by no more than a factor of two, 
the errors in the measurement of dip angle would be halved, and the energy 
resolution considerably improved. 


Thus, in the fields referred to above, as in many other fields of investiga- 
tions employing the nuclear emulsion technique, the desirability of inhibiting 
shrinkage is now well recognised. The advantages of shrinkage inhibition 
could be listed as follows: 


1. Angles of dip of tracks can be measured with accuracy comparable 
to that with which the horizontal angles in the plane of the emulsion can be 
measured. 


2. Vertical projections of tracks are more than in shrunk emulsions, 
and hence range measurements can be made with improved accuracy. 


3. As compared with a shrunk emulsion, the average background grain 
density is reduced by a factor which is the ratio of the shrinkage factors of 
the shrunk and shrinkage inhibited emulsions. This reduction in background 
will be of particular value in the range measurement of low energy tracks— 
protons around 0:5 MeV and low energy electrons. 


4. Accidental coincidences between distinct tracks is minimised. 
PREVIOUS METHODS 


After the washing stage, and prior to drying, it is customary to soak 
nuclear emulsions in a 5% solution of glycerine, to minimise strains 
which are set up in the emulsions during the drying process. It was 
found by Meulemans and Mignone* that a concentration of 10-11% of 
glycerine gave a shrinkage factor close to unity. Glycerine however is not 
the best shrinkage inhibiting agent because (1) its refractive index is slightly 
different from that of gelatine and the immersion oils, and hence the optical 
quality of the emulsions is impaired, (2) it is hygroscopic—the glycerine treated 
AS 
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plates must therefore be stored and analysed at constant relative humidity; 
this is a major shortcoming of this technique, and (3) the emulsion becomes 
soft and tends to melt under the intense beam of light in the microscope. In 
spite of these difficulties, this method has been employed, for the sake of 
improved energy resolution, in the study of the fission neutron spectrum of 
U2 by Cranberg et al.4 


Another method has been developed by Dansyz and Yekutieli,> in which 
after the washing stage, the water in the emulsion is replaced by alcohol, 
by immersion in a succession of baths of gradually increasing concentration 
of alcohol. The plates are then soaked in a solution of 30 gm. of resin in 
100 c.c. of alcohol. A 300 micron emulsion thus soaked for a period of 
48 hours gave a shrinkage factor close to unity. 


Our METHOD 


A systematic examination of the various substances available for 
incorporation into the emulsion revealed that canada balsam was by 
far the best substance. It has a refractive index same as that of the 
emulsion, and unlike glycerine it is found to insulate the emulsion from 
humidity variations of the surroundings. It is only partially soluble in alco- 
hol, and hence cannot be incorporated in the emulsion by the method of 
Dansyz and Yekutieli. It is however very soluble in benzene which in turn 
is miscible with alcohol, but not water. We did a preliminary experiment 
to ensure that benzene did not damage a processed emulsion. The following 
procedure was then adopted to impregnate the emulsion with canada balsam. 


After washing, the plates were leached of water by immersion in success- 
ively increasing concentrations of alcohol in water. The alcohol concentra- 
tion was increased in 10% steps, and the plates were maintained in each stage 
for a period of time which depended upon their thickness—75 minutes for 
200 micron plates and 6 hours for 400 micron plates. After the 100% alcohol 
stage, the plates were immediately put in a 10% solution of benzene in alcohol. 
They were then taken through benzene-alcohol mixtures, with the benzene 
concentration increasing in 10% steps, the time for each stage being the same 
as that for the water alcohol mixtures. After the 100% benzene stage, the 
plates were immediately transferred into a solution of canada balsam in 
benzene, the optimum concentration of which was arrived at after a series 
of trials. A concentration of 100 gm. of canada balsam in 25 c.c. of benzene 
with a soaking time of 12 hours for 200 micron emulsions (Ilford C2), gave 
a shrinkage factor close to or lesser than unity. After removal from the 
canada balsam, the plates are allowed to dry and then cleaned with benzene 
soaked in cotton wool. The actual value of the shrinkage factor depends 
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not only on the concentration of the canada balsam-benzene solution, but 
also on the type, history and pretreatment of the plates. This aspect will 
be discussed separately. The times for the water-alcohol and alcohol-ben- 
zene stages, and the soak in the canada balsam-benzene solutions are indi- 
cated in Table I. 





TABLE I 
Thickness in microns a 100 200 400 600 1000 
Water-alcohol or alcohol- - 
benzene stage .. .. 35min. 75 min. 6 hrs. 10 hrs. 12 hrs. 
Benzene-canada balsam 
solution wv .. 6hrs. 12 hrs. 60 hrs. 100 hrs. 120 hrs. 





The temperatures of all the above solutions are kept below 10°C. In all 
stages of the processing and particularly in the canada balsam stage it is 
important to maintain the plates in well levelled trays. 


DETERMINATION OF THE SHRINKAGE FACTOR 


To establish the reliability of this method of shrinkage inhibition, a 
simple yet accurate method has been employed for the determination of the 
shrinkage factor. Before processing, the thickness of the emulsion along with 
the glass backing is measured at six different points on the plate, using a dial 
gauge. The measurements are accurate to within one micron. The thickness 
of the processed emulsion with the glass backing is again measured at the 
same points. The thickness of the processed emulsion only at these points 
is then measured with a microscope. From the set of three measurements 
at each point, the thickness of the unprocessed emulsion and hence the 
shrinkage factor at that point can be obtained. 


For any particular plate, the value of the shrinkage factor at six different 
points is thus obtained, and from these, the mean value of the shrinkage 
factor for the whole plate is obtained. This mean value which lies around 
unity for the shrinkage inhibited plate, has a standard deviation of 0-03. In 
the case of some plates which were not held horizontal during processing, in 
particular during the canada balsam stage, the standard deviation for the 
shrinkage factor of a plate was as high as 0:07. Experiments have however 
shown that by carefully controlling conditions, it is possible to achieve uniform 
impregnation of canada balsam, and thus obtain values of shrinkage factor 
with standard deviation as small as 0-01. In a typical experiment six Ilford 
200 micron C2 plates belonging to the same batch and hence with the same 
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history, were identically processed and subjected to the canada balsam treat- 
ment. The thickness of the unprocessed emulsions varied from 187 microns 
to 228 microns. The variations on individual plates were much less—189 to 
204 and 220 to 228 in two typical cases. The values of shrinkage factors 
for individual plates, obtained as a mean of six independent determinations 
at six different points on each plate, were 0-80 + 0-02, 0-75 + 0-01, 0-78 
+ 0:01, 0-81 + 0-02, 0-79 + 0-01 and 0-78 + 0-02. The average value 
of the shrinkage factor for all the six plates is 0-79 + 0-02. Four other 
Ilford C2—200 micron plates which were processed and dried in the conven- 
tional manner (alcohol drying without glycerine) gave shrinkage factors of 
2°53 + 0:03, 2:51 + 0-02, 2:-49+ 0-03, and 2-61 +0-03. The average 
value for the four plates was 2:52 + 0-06. It will thus be seen that a plate 
treated with canada balsam yields a value of shrinkage factor which is con- 
stant within experimental errors over the whole of the plate. 


As indicated earlier, the value of the shrinkage factor differs from plate 
to plate, depending upon the type, history and pretreatment. For example, 
various 200 micron plates processed and shrinkage inhibited by identical 
processes, gave shrinkage factors varying from 0-7 to 1-2. The important 
factor however, was the constancy of the value of the shrinkage factor over 
the whole of any individual plate. Therefore, in any shrinkage inhibited 
plate, it will be sufficient to make the various thickness measurements at say, 
three points. The emulsion thickness before processing at these points is 
thus known, and the shrinkage factor of the emulsion at any time when track 
measurements are being made, can be determined by measuring the emulsion 
thickness at these points under the microscope. We can thereby eliminate 
errors due to minor variations in the thickness of a plate due to various factors 
like time, storage conditions, etc. 


In order that gamma ray measurements with heavy water loaded emul- 
sions may be improved further, an attempt was made to increase the final 
thickness of the emulsions by increasing the concentration of canada balsam. 
Four plates from the same batch were identically processed and shrinkage 
inhibited in the following manner. Two were placed in a solution of the 
already recommended concentration, viz., 100 gm. of canada balsam in 25 c.c. 
of benzene. The other two were placed in a solution of 150 gm. of canada 
balsam in 25 c.c. of benzene. The former two plates gave a shrinkage factor 
of 1:07 + 0-02, while the latter two gave 1-08 + 0-03. It will thus be seen 
that any sizeable reduction in the value of the shrinkage factor cannot be 
achieved by increasing the concentration of canada balsam. The exact 
nature of the dependence of the value of the shrinkage factor on the history 
and pretreatment of the plates has not been investigated. 
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UNIFORMITY OF SHRINKAGE INHIBITION IN DEPTH 


To ensure further, that the shrinkage inhibition was uniform throughout 
the depth of the emulsion the following investigation was undertaken. The 
ratio of the number of background grains in equal areas of the shrunk and 
shrinkage inhibited emulsions (from the same batch and processed identically) 
was determined at various depths. If the shrinkage inhibition is uniform, 
this ratio should be the same for all depths, and should be equal to the ratio 
of the shrinkage factors of the shrunk and shrinkage inhibited emulsions. 
In a typical study, the ratio of the shrinkage factors of 2 plates was 3-11. 
The ratio of the grain counts made at three different depths gave the values 
3-46, 3-28 and 3-29. In addition, emulsions soaked in benzene-canada 
balsam solution for times much larger than those given in Table I, e.g., 18 
hours for 200» emulsions and 90 hours for 400 » emulsions, did not have 
shrinkage factors different from those emulsions which had been normally 
treated. 


VARIATION OF THICKNESS WITH HUMIDITY 


In view of the hygroscopic nature of glycerine, it would be expected that 
the thickness of a glycerine impregnated emulsion would be sensitive to 
humidity variations of the surroundings. A canada balsam treated plate, 
however, should be expected to render the plate almost immune to such 
humidity variations. In order to study the behaviour of differently treated 
emulsions under conditions of varying humidity, the following experiment 
was undertaken. 


Five plates were treated as follows: The first plate was removed from 
the 100% alcohol stage and allowed to dry. The second plate after removal 
from the alcohol stage, was subjected to a 5% glycerine treatment. The third 
plate was treated similarly to plate 2, except that 15% glycerine was used. 
Plate 4 was removed from the 100% benzene stage and allowed to dry. Plate 5 
was subjected to the canada balsam treatment. All the five plates were 
then placed in a humidity chamber. The humidity in the chamber was held 
constant at any particular value by using solutions of glycerine in water of 
suitable concentrations.® At the various values of humidity, emulsion: thick- 
nesses were measured under the microscope, at three well defined points on 
each plate. Care was taken to ensure that the emulsion thickness at any 
particular value of humidity had reached equilibrium. This was done by 
maintaining the emulsions at a particular humidity value until thickness 
readings taken at any time were reproducible atleast 48 hours later. These 
readings were then taken as the thicknesses corresponding to the particular 
value of humidity, . 
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The results obtained are represented in Fig. 1. It will be seen from the 
graph that whereas the glycerine soaked plates are extremely sensitive to 
humidity variations, the other plates are not. Table II gives the ratios of 
the thicknesses at the two extreme values of humidity at which measurements 
were made. 
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fact that gelatine is itself hygroscopic and hence the addition of canada bal- 
sam improves its thickness preserving qualities by inhibiting water absorp- 
tion. It was further observed that the optical properties of glycerine treated 
plates deteriorated progressively with increase in humidity. At values of 
humidity above 70% the glycerine plates were misty especially in the lower 
regions, and observations under the microscope were rendered extremely 
difficult. At low humidities, the alcohol and benzene dried plates tended 
to peel off from the glass backing, whereas the glycerine and canada balsam 
treated plates remained completely unaffected. 


DISCUSSION 


From the above findings, it is obvious that the canada balsam treated 
plate has distinct advantages over the glycerine treated one. The canada 
balsam treated plate can be stored and measured in a laboratory where faci- 
lities do not exist for controlling humidity. The uniformity of the shrinkage 
factor value over the whole of a canada balsam treated plate having been 
established, it should now be possible to make accurate measurements to 
determine the mean angle of scattering of a track by vertical deflections, by 
the method of Mabboux.? Further, it should be possible to work with rela- 
tive ease with emulsions in which foreign bodies like wires and capillaries 
have been incorporated. Stripped emulsions processed without glass backing 
could be shrinkage inhibited with canada balsam and the emulsion could 
then be rolled on to a glass plate, thus ensuring perfect optical contact between 
emulsion and the backing glass. 
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A RELATIVISTIC irreducible wave equation for a particle having the two mass 
and spin states (m, 3/2) and (Am, 4) and positive-definite charge was proposed 
by Bhebha.! It can be rewritten? in a form analogous to that of the Rarita- 
Schwinger*® equations for particles of half odd-integral spins. The purpose 
of this paper is to write down the non-relativistic approximation of the equa- 
tions in an external magnetic field, with a view to determining the magnetic 
moment of the particle in the two states. It is convenient to use the Rarita- 
Schwinger form of the theory for this purpose, since only the Dirac matrices 
occur in the equations in this form, rather than the complicated 20x20 
matrices of the original theory. We follow the method of Tiomno* who 
obtained the non-relativistic approximation of the Fierz-Pauli spin 3/2 equa- 


tions. The magnetic moment in the spin 3/2 state is found to be e/3mz, 
which is also the value for the pure spin 3/2 case. The value for the spin 
1/2 state is e/2Am {1 + 2Ad?(1 + A)}, and differs from that of a Dirac 
particle. 


1. The equations of motion in the presence of an electro-magnetic 
field are 


(Pry? + m) dy — 4(Prvk + Prva) o! + 4 7 (Pry” — m) yo¥® + dPpp = 0, (1) 


(Pry” + Am) + — Ppp" = 0, (2) 
where 


Pi) 
P oo i — é. 
b= i Ax, 


and A; is the 4-potential of the external field. Multiplying (1) by y* and 
Pk and summing over k we obtain, respectively, 


% Pr? — 4 myzpf" + dPry ys = 0 (3) 
and 


= Pry Py — 3 Pry ys + mPph + ie Fimy's™ 


— tie Fumy'y™yr" + dP% = 0, (4) 
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where 


— dsAm dA] 


‘oxt — axm 
Multiplying (3) by Psy* and subtracting from (4), we get 
m Pr? = — ie Fumy'b™ + 3 ie Fum y'y™yr" — Cie Fimy'y™. (5) 


From (2) and (3) we have 


] 
Po” = (95 aqay (mvei” + 3ddomp3 6) 
and 
| 
aa " 
Equations (5) and (6) now give 


(2 + 3d) , 
m 


yl me —~ Sl 4 {— ie Fimy'b™ + 2 ie Fun y'y™y 


ad. 
— 5 ie Fim yy}. (8) 


2. In the non-relativistic approximation the equations of motion should 
involve only two independent components for the spin 1/2 state and four inde- 
pendent components for the spin 3/2 state. We shall write down the equa- 
tions of motion for the two cases as Schrédinger equations, the Hamiltonian 
being correct to terms of order 1/m. One can then remove the small compo- 
nents by means of a Foldy-Wouthuysen transformation. 


From (8) we have, to terms of order (1/m)?, 


8 
yr? = — 3dnp + OTIS I e_ Braytym — 201 +2) Frmyy™. (9) 


We now wish to calculate Py (y,") to terms of order 1/m. For this purpose 
we assume that the external field is a constant, weak magnetic field, 
H = ( — Fis, — Fes, — F;) and that 
[Prt] << my}. 
Note also that, substituting for Pt” from (6), equation (1) can be written as 
ke A+) kyr Am _ ky a pky yt 
(Pry" + m) (2 + 3d?) my" y rip (2 + 3d) yp +P Yr 


+ 4 y*Ppy" (yep5) + dPhy = 0, (10) 
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Making use of (7), (10) and (9) we obtain, to terms of order 1/m, 
Pay? (yet) = 3dntmp + CEA — 34° — 2) ayy 


a 3iedA 


mF imy'y™b. (11) 


Substituting for y,%" from (9) and for P,;y" (ys#5) from (11), equation (1() for 
k = 1, 2, 3 becomes 


(Pyy™ + m) bk + ddm (1 + dA) phy oa i (1 + 3d?A) YF imylyi™ 


+ om (1 + d? — A+ Ad®) YF ymy'y™p = 0. (12) 


Similarly, making use of (6) and (9), equation (2) becomes 
_ ied , ted? (1 + A) 
(Pry” + Am) $ + 7 Fimy'p™ + et = Fimy'y™ = 0. (13) 
3. Let us now consider the spin 3/2 state. In the force free case, in 


this state % is identically zero. In the presence of an external field we must 
therefore have 


b= a Frmy'$™, (14) 


where a is a constant. Substituting this expression for % in (13) and making 
use of (12) one easily obtains 


sae ied | 
m* (1 + A) 
Equation (12) now becomes 
(Pry™ + m) oF + i Y*Fimy'y™ = 0, (15) 


neglecting terms involving squares of the field strengths. This is just the non- 
relativistic approximation of the pure spin 3/2 equations obtained by Tiomno, 
and has been further reduced by him, by means of a Foldy-Wouthuysen trans- 
formation, to the form 


dg 


: > 
l 


at 1 he A)? 4 A e > 
dt 13m (P — ¢ , Te o— 3,2 - Hf 4, 








Ww = 


Mise 
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where ¢ is a four component quantity. The three components of z are the 


3x3 matrices representing an angular momentum of 3/2. The magnetic 


moment of the particle in the spin 3/2 state is thus e/3m =. 


4. For the spin 1/2 case the wave-functions %* can be determined from ws. 
In the force free case we have, in the non-relativistic approximation [see 
equation (21 a) of Aj, 


pe — — dryhy, (17) 
In the presence of an external field this relation will be of the form 
gk — — dryk + terms involving Fim. (18) 


Substituting (18) in (13) we obtain 
(Pry” + Am) % + vo 7 — a) Fimy'y™p = 0. (19) 


This may be rewritten as 


) a? 
j™ = for, + (p— eA) yy — dmy— 529 (1 — 2) imply 4. (20) 


By means of a Foldy-Wouthuysen transformation this equation can be trans- 
formed to 


.2 o > > 
j>Y — {— rm + edo — 4K (p — eA)? 


ya (PUR9 SI. an 


For the positive energy states only the first two components of ¢% are non-zero. 
Making the tranformation 


b — eidmty, 
we obtain, for positive energy states, 
; i a ae 2 ge == a eo +> 
is? = {eho + ayq(P — eA)* — H- ( + sy, )ohy, (22) 


being now a two component spinor. The particle thus has the magnetic 
moment 


+ 
M = 


e 2 per 
2m {1 + 2Ad? (1 — A)}o, (23) 
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Note that it reduces to the value for a Dirac particle if A = 1, i.e., if the masses 
in the two states are equal. 


It is a pleasure to thank Dr. M. K. Sundaresan for stimulating discussions. 
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1. INTRODUCTION 


THE study of the daily variation of cosmic ray intensity has grown in import- 
ance in recent years with increasing realisation that it provides a powerful 
tool for understanding the anisotropy of primary cosmic radiation as well 
as the electromagnetic state in interplanetary space. With the solar daily 
variation of cosmic ray intensity which arises due to the spin of the earth, it 
is possible to investigate the anisotropy in various directions by making 
suitable directional studies at different latitudes. A mass of data collected 
during the last three or four years as well as a re-examination of earlier data 
have brought to light numerous experimental facts which have been reviewed 
recently by Sarabhai et al.1 The relevant facts may be stated as follows :— 


1. The daily variation cannot at all times be satisfactorily described 
in terms of the diurnal component alone. It has been observed, however, 
that in general the first two harmonic components are adequate to represent 
fairly satisfactorily the main features of the daily variation. 


2. The 12-monthly (annual) mean daily variation of meson intensity 
exhibits large long-term changes of amplitude and of the time of maximum 
of the diurnal as well as the semi-diurnal components.? 


3. During several years the 12-monthly mean daily variation, parti- 
cularly at low latitudes, exhibits two maxima instead of one (Fig. 8 of Ref. 2). 
Sarabhai et al.? have tried to show that the 12-monthly mean daily variation 
can be explained in terms of two types of variations, one having a maximum 
during the day at about 1100 to 1300 hours. and the other at night at about 
0100 to 0300 hours. The relative contribution of these two types of varia- 
tions changes from year to year and the data collected over the past twenty 
years is suggestive of a relationship of these changes with the 22-year cycle 
of solar activity. This makes the maximum of the diurnal component of 
the daily variation shift by as much as 10 to 12 hours over a number of years. 
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4. When data on individual days are examined, there is a tendency for 
the time of maximum of the diurnal component to occur either at noon or 
during the early morning hours.’ Such days usually occur in groups that 
sometimes recur after 27 days. Some of the groups appear to be associated 
with active solar regions. 


5. High values of Cp, daily geomagnetic planetary index, are normally 
associated with days having large amplitudes.’ However large amplitudes 
are not always associated with magnetically disturbed days or the central 
meridian passage of solar regions of visible activity.® 


The above facts lead to the conclusion that the daily variation and the 
anisotropy of cosmic radiation either on individual days or averaged over a 
group of days is of a highly variable character. Our knowledge of the per- 
manent anisotropy, if indeed it exists, is meagre at the present moment. 
Hence for the interpretation of the anisotropy it is appropriate for the time 
being to take into account only the variable anisotropy about which we have 
a number of well established facts regarding its form as well as its solar and 
terrestrial relationships. 


In the past few years, several attempts have been made to explain some 
of the average characteristics of the diurnal component of the daily variation 
on a model involving the modulation of primary cosmic ray intensity by 
changes in the electromagnetic condition in interplanetary space. Brunberg 
and Dattner? have tried to explain an anisotropy which would produce a 
daily variation of average amplitude of 0-3% and a time of maximum near 
noon. They suggest that a tangential anisotropy in the 18-hour direction 
might be created in the solar system by the rotation of the magnetic field of 
the sun. Alfven® has moreover shown that the disturbance in this field by 
the solar emission of ionised matter can give rise to a radial anisotropy 
which would alter with solar activity and could therefore account for the 
11 years variation in the time of maximum of the diurnal component shown 
by Sarabhai and Kane’ and Thambyahpillai and Elliot.” 


To explain the anisotropy on days when there is a magnetic storm, 
Nagashima" has studied in detail the implications of the suggestion made 
by Elliot and Dolbear'* and Alfven?* that the storm time anisotropy is due 
to the presence of a polarised beam of ionised matter in the vicinity of the 
earth. His analysis shows that this anisotropy in the primary radiation is 
directed along the 12-hour direction and explains for the first time the increase 
in amplitude and advancement of the phase of the diurnal variation on mag- 
netically disturbed days. 
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Apart from the difficulty that arises due to the daily variation exhibiting 
characteristics markedly different from those assumed by Brunberg and 
Dattner,? none of the theories advanced so far accommodate the existence 
of the two main types of variable anisotropies which are responsible for the 
day and night maxima in the daily variation of cosmic ray intensity at 
low latitudes. In this paper an attempt has been made to understand the 
variable character of anisotropy by modulation of the primary cosmic ray 
intensity by ionized solar-matter as it is ejected from the sun carrying a frozen 
magnetic field. The anisotropy is derived from the measurements of the daily 
variation of meson intensity at low latitudes. 


Il. CHARACTERISTICS OF THE SOLAR CORPUSCULAR EMISSION 


The recurring disturbances in the earth’s magnetic field have been attri- 
buted by a number of workers" to the presence of streams of ionized matter 
emitted by the sun. Although a 27 days recurrence of these disturbances 
over two to three solar rotations suggests that the emission is active over a 
long period, individual disturbances last for a few days. This indicates that 
the earth is engulfed in the streams for only a few days. It is, therefore, 
assumed that a mass of gas ejected by the sun may be looked upon as a 
beam with parallel sides perpendicular to the plane of the ecliptic. 


Alfven!® has shown that a highly conducting medium moving out of a 
region having a high magnetic field tries to drag the field with it. A mass 
of gas which is expelled by the sun is ionized and he has suggested that it will 
carry with it the solar magnetic field. Recent work on the general magnetic 
field of the sun by Babcock and Babcock’* puts a limit of | gauss to its mag- 
nitude and suggests that it is concentrated near the polar regions. In solar 
equatorial latitudes there are intense local magnetic fields associated with 
active regions on the sun. It is reasonable to expect that beams along the 
plane of the ecliptic are formed by solar corpuscular emission from solar 
equatorial latitudes and they would carry therefore a trapped field which is 
not derived from the general solar dipole field but from the local magnetic 
fields in the neighbourhood of active regions from which the ionized matter 
is supposed to be emitted. On this basis it would be possible to have the 
trapped field having no preferential orientation with respect to the solar 
dipole field. For simplifying calculations it is assumed that this field is uni- 
form within the beam. Swann?’ has explained how such uniformly magnet- 
ised beams with no external magnetic field can be realised in practice. 


The trapped magnetic field moves with the beam with the same velocity 
as the solar corpuscles and in a co-ordinate system moving with the cor- 
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puscles, there is only a magnetic field present. But an observer on the earth, 
with respect to which the magnetic field is in motion, observes in addition 
to a magnetic field a charge displacement in a direction perpendicular to the 
velocity of the beam and the magnetic field. The beam therefore appears 
to have crossed electromagnetic fields and cosmic ray particles traversing the 
beam gain or lose energy depending upon their orientation with respect to 
the electric field. An isotropic distribution of primary particles is dis- 
turbed and an anisotropy is created in some directions in space. 


III. ANISOTROPY OF THE PRIMARY COSMIC RADIATION* 


In evaluating the effect of a beam of solar ionized matter on the aniso- 
tropy of primary cosmic radiation we neglect the inclination of the axes of 
the sun and the earth with the ecliptic. For observing stations at low lati- 
tudes, the primary cosmic ray particles incident from the vertical come from 
directions close to the equatorial plane, and hence only the acceleration or 
deceleration of particles crossing the beam in this plane is of consequence. 
This is produced by an electric field which develops due to the component 
of the trapped magnetic field in a direction perpendicular to the plane of the 
ecliptic. With different orientations of the magnetic field within the beam, 
this component could be reversed in direction. In the present work, the 
effect Of this component alone on the anisotropy of cosmic radiation is in- 
vestigated, and the results are most closely applicable to observations at the 
equator with narrow angle telescopes pointing to the vertical. However 
omnidirectional instruments, such as the ionization chamber, measures in 
addition primary particles inclined to the equatorial plane. In this case, 
as well as for observations at high latitudes, the contribution of the compo- 
nent of magnetic field parallel to the ecliptic would also have to be considered. 


Fig. | shows the cross-section of the beam in the plane of the ecliptic. 
The ionized matter in the beam is moving with a velocity v along the Y-axis. 
The beam has a width equal to 5 along the X-axis. It is supposed that the 
beam is situated symmetrically with respect to the plane of the ecliptic. The 
magnetic field trapped within it has a component H,’ along the Z-axis and is 
called positive along the positive direction of the Z-axis. To a stationary 
observer on the earth the beam will have a component of magnetic field 
Hz = H,’/(1 — v?/c?)t along the Z-axis and an electric field (v/c) Hz along the 
X-axis. The potential difference between the two parallel sides of the beam 
is therefore (v/c) bHz and cosmic ray particles of charge e crossing such a 
beam will change their energy by (v/c) ebH3. 


* Part of the analysis in this section is reproduced from the original paper by Nagashima." 
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Fic. 1. Cross-section of the beam in the plane of the ecliptic. 


The direction of the velocity of a particle in space can be uniquely deter- 
mined by two angles ¢ and # where ¢ is the inclination to the plane of the 
ecliptic (XOY plane in Fig. 1) and y# is the inclination in the plane of the 
ecliptic to the sun earth line (Y-axis in Fig. 1). If Wj, (¢,, %,) represent the 
energy and the direction of the particle before entering the beam and W,, 
(¢., ,), those of the particle after crossing the beam we get : 


Since the velocity of ionized matter in the beam is much less than the velocity 
of light, v<c. Also if 4 Wy < W, the direction of the particle with respect 
to the plane of the ecliptic is not appreciably altered. In consequence 
db, © dy 

The particle is however deflected in the plane of the ecliptic and as shown in 
Fig. 1 if ¢ is not very large 

COS p. = Cos #, — D/py (2) 

= cos %, — sec ¢ (ebHz)/W, 


where py is the projection of the radius of curvature on the plane of the 
ecliptic. 

It is seen from Eq. 2 that if b > 2py, cos 4. << — 1 and #, is imaginery. 
Such particles are not able to cross the beam but after entering the beam are 
deflected so as to emerge with unchanged energy on the same side of it from 
which they entered. Particles having b < 2p, can cross the beam only if 
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they enter it in favourable directions. The energy W of a proton which can 
cross the beam of width 6 in a direction ¢ is given by 


W = W,,,, (¢, 6) = 4 ebH; sec 4 (3) 

In Fig. 2 have been shown the deflections in the plane of the ecliptic of 
particles having energy less than W,,,, (0, 5). It is noticed that if initially 
there is an isotropic distribution of such particles, the deflections within the 
beam do not alter the intensity of these particles in any direction outside the 


beam. Thus the isotropy of cosmic radiation of energy less than W,,,, (0, 5) 
is not disturbed by the beam when an observer is outside it. If, on the other 
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Fic. 2. Deflections of low energy cosmic ray particles W < Wmin(0, 5) in the magnetic field 
within the beam. 


hand, the earth is within the beam, the particles have to cross only a part of 
the beam to reach the earth. Since the strength of the magnetic field is the 
same as before but the width is reduced, the minimum energy of particles 
which can reach the earth is reduced and the anisotropy therefore extends 
to particles with energy less than W,,,, (0, 5). When the earth is at a distance 
x from one side of the beam, the minimum energy of the anisotropic radiation 
Win (0, x) is given by 


Win (0, x) = 4 exH;z or 4.e (6b — x) Hz (4) 
whichever is smaller. 


The anisotropy of cosmic rays also depends on the position of the earth 
with respect to the beam. The distance between the earth and the beam 
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changes due to the spinning of the sun and the motion of the earth along the 
ecliptic. Beams approach the earth in the same direction as the motion of 
the earth, but soon overtake it. We have to consider three cases; firstly 
when the beam approaches the earth, secondly when it envelops the earth 
and finally when it recedes from the earth. In Fig. 3 we show the effects in 
each case when H, is positive. When the observer is on the right side of the 
beam, particles of energy W after crossing the beam from left to right gain 
an energy 4 Wy, and have their final velocity directions lying in the shaded 
region in Fig. 3(a@). This region of acceleration extends from %p** to 7 
where #p** is given by Eq. 2 as 


acc —— 1 — 2 sec @® —ain 
bp cos~ (1 2 sec > WwW (5) 
There will be an increase of intensity of particles arriving in directions lying 
in this region, while those arriving from other directions either have not passed 
through the beam or, if passed through, have not crossed the beam and 
therefore do not undergo any change in their energy or intensity. 


In Fig. 3 (5), the earth is enveloped in the beam but is at a distance x 
from the left side of the beam. The beam can now be considered as made up 
of two parts, one having width equal to x and the other (b — x). Particles 
arriving from the left side are accelerated and gain an energy by 
4 W. = (v/c) exH, and those coming from the right are decelerated and lose 
energy by 4 Wip_x) = (v/c) e(b — x) Hz. Due to the deflection of parti- 
cles within the beam, the region of acceleration extends from ¥, to pz) 
while that of deceleration from ¥%,p_y, to %, where 


2 = cost (1 — 2sec ¢ pW) and 7 > bz >0 (6) 


tib-x) = cos (- 1+ 2sec¢ b + x Ma) and 27 > do-2) > 7 (7) 





When the earth is in the centre of the beam, i.e., x = 5/2, Yy and pz) 
are separated by 180°. We thus see that when the earth is within the beam, 
all directions in space, as viewed from the earth, are divided into either accele- 
rated or decelerated regions. The intensity of the particles will therefore 
be more or less according as the instrument points to the region of accelera- 
tion or deceleration respectively. 


In Fig. 3 (c) the earth is on the left side of the beam and particles of 
energy W while crossing the beam lose their energy by 4 Wy, and appear in 
a shaded region bounded by #p** and 27 where 
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doit = cost(— 1+ 2secp an) and 2e >We >2 


In this region the intensity of the particles is below normal. 


CONDITIONS FOR POSITION OF THE EARTH 
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Fic. 3. Deflections of primary particles of energy 2W»min (0, 5) and anisotropic region when 
the earth is (a) on the right side, (6) within and (c) on the left side of the beam. The component 
of the trapped magnetic field H, is positive and the type of daily variation as a function of local 
time is measured by an instrument on the spinning earth in the absence of the geomagnetic field. 

Because of the spin of the earth, a fixed telescope at a low latitude station 
would scan all directions in the plane of the ecliptic and would thus view 
during each day the regions where the isotropy of cosmic radiation is dis- 
turbed. The intensity measured by the telescope would therefore exhibit 
a solar daily variation. Taking the walls of the beam along the sun-earth 
line, i.e., 12-hour direction and neglecting the influence of the earth’s mag- 
netic field, we show in Fig. 3 4 W, as a function of local time T cn which is 
dependent the type of daily variation to be expected in each of the three cases. 
There is an increase of intensity from 12 hours to (27 —. #p*)/15 hours in 
Fig. 3 (a’), while a decrease from 00 hours to (27 — pp**)/15 hours in 
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Fig. 3 (c’). The time of occurrence of the increase in case 3(a’) and the 
decrease in case 3(c’) is separated by 12 hours. The maximum of the daily 
variation therefore occurs between 12 and 18 hours irrespective of the relative 
position of the beam and the earth. When the earth is within the beam, 
this maximum, as shown by Nagashima," is shifted towards 12 hours but is 
still between 12 and 18 hours. 


We can summarise the influence of a beam with a positive component 
of the trapped magnetic field as follows. From the time when the beam is 
effective while approaching the earth, a region in space in which particles 
are accelerated is formed and this results in an increase in mean intensity in 
the afternoon. As the earth enters the beam, a new region of deceleration 
is created. This region goes on progressively decreasing but with increase 
in the amount of deceleration produced. The region of acceleration simul- 
taneously increases but with a decrease in the amount of acceleration. The 
mean intensity therefore progressively decreases and attains its original value 
at the centre of the beam when both the regions are equal in width. As the 
earth leaves the beam, the region of acceleration disappears leaving behind 
a decelerated region and a decrease of mean intensity which returns to its 
normal value when the beam ceases to be effective. 


In Fig. 4 we illustrate all these cases when H, is negative, i.e., directed 
opposite to what has been considered so far. The electric field is now directed 
from right to left so that in Fig. 4 (a) when the beam is approaching the earth 
we get a decelerated region from 0 to #p%** where 


pp = cos (- 1 + 2sec¢ Wan ?)) and 7 > pp > 0 (9) 


In Fig. 4 (5) there is an accelerated region from #z to %p_~, and a decelerated 
region from ¥,p_z) to #z where 


x W,,,, (0, b 
px = cos} (- 1+ 2sec¢ b San) and 7>7%,>0 (10) 
and 
b— xW,,, (0, 6 
%(b-x) = cos (i- 2 sec a : Mo) and 27 > $o-x >a (11) 


In Fig. 4 (c) when the earth leaves the beam we get an accelerated region from 
7 to %p** where 


pp = cos! (: — 2sec d wo) and 27 > p** >T (12) 
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Fic. 4. Deflections of primary particles of energy © 2Wmin (0, 5) and anisotropic region 
when the earth is (a) on the right side, (6) within and (c) on the left side of the beam. The com- 
ponent of the trapped magnetic field H, is negative and the type of daily variation as a function 
of lecal time is measured by an instrument on the spinning earth in the absence of the geomagnetic 
field. 

The type of daily variation recorded by an instrument on the spinning 
earth, neglecting the effect of the geomagnetic field, is illustrated in 
Fig. 4 (a’, b’, c’). When the earth is approaching the beam there is a decrease 
of intensity from (27 — Yp**)/15 hours to00 hours and when it leaves the 
beam there is an increase of intensity from (27 — p**)/15 hours to 12 hours. 
The maximum of daily variation, when various energies of primary radiation 
are considered, is therefore expected to be between 06 and 12 hours. 


The electromagnetic fields of the solar beams thus create anisotropy in 
the primary radiation which when the trapped magnetic field has a positive 
component gives rise to a maximum in the daily variation between 12 and 
18 hours. With a reversed component of the magnetic field this maximum 
occurs between 06 and 12 hours. The times of maxima relate to the condi- 
tion when the bending of the trajectories of cosmic ray particles in the geo- 
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magnetic field has not been allowed for. In the former case, when the beam 
is effective we should get an increase of mean intensity followed by a progres- 
sive decrease while in the latter case there should be a decrease of intensity 
followed by an increase. 


IV. DrURNAL VARIATION OF MESON INTENSITY PRODUCED BY THE 
ANISOTROPY OF PRIMARY COSMIC RADIATION 


The foregoing analysis gives the nature of the anisotropy that will be 
created by the electromagnetic fields in solar beams. The observable effect 
of the anisotropy in terms of a solar daily variation is however measured on 
cosmic ray secondaries produced by the primaries after they have suffered 
deflection in the geomagnetic field and interacted with matter in the atmo- 
sphere. Since the anisotropy, the deflection in the geomagnetic field as well 
as the interactions in the atmosphere are dependent on the energy of the 
primary particles, it is preferable for the purpose of evaluating the daily 
variation to consider the entire cosmic ray spectrum above the geomagnetic 
cut off rather than a mean energy of the primary radiation, as is usually done. 


When the primary radiation is isotropic the intensity of cosmic radiation 
I(A, 4) observed at an atmospheric depth h g./cm.? and at geomagnetic lati- 
tude A is given by 

=27 @=7/ co 
I(,h)= ff J m(E,h, 0,0)i(E)F(6,0)sin@dad@dE (13) 
a=0 6=0 Ej, 0,a 
F (8, a) denotes the sensitive area of a recording apparatus at a zenith angle 
§ and azimuth angle a and E,, 4, , is the corresponding cut-off energy at the 
geomagnetic latitude A. m(E,h, 8, a) and i(E) dE are the overall multiplicity 
and differential energy spectrum respectively of the primary radiation. 


If by the presence of the beams, primary particles of energy E gain an 
energy 4E = 4W for a certain period of time, the change in the intensity 4I 
is given by 

a=2m Q=T/2 


a=7 3 


J mE, h, 0, a) iE) g (B) AE (Tz, 9, «) F (0, a) 
E 


x sin 8 da dédE (14) 


where Cy = Ey, 9, OF Wain (¢, x) whichever is greater. The factor g (E) 
is introduced to take into account the change in the differential energy spec- 
trum of primary protons due to the additional energy gain and is given by 
Fonger?® as 


g(E) = _ Fesieaoe 4 B| where 8 = 2-07 (15) 
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For meson intensity in the lower atmosphere, when h > 700 g./cm.*, 
m(E, h, 8, a) can be approximated as m(E, h)cos?6. The local time Ty, , , 
refers to the direction of the particle outside the earth’s magnetic field. Its 
relation with the local time fg, 9, . on the surface of the earth is given by 


Tr, 9, a = te, 9,a + ¥z, 9,0 (16) 


where #9, is an angle through which a cosmic ray particle of energy E 
arriving at a zenith angle @ and azimuth angle a is deflected in the east-west 
plane (the plane of the ecliptic in the present work) in the earth’s magnetic 
field. These values can be read out from the curves supplied by Brunberg 
and Dattner.1® To evaluate the change of intensity at the surface of the 
earth Ty, 9, q in Eq. 14 is replaced by tg, 9,4, To obtain the diurnal varia- 
tion, the values of AI/I are calculated by numerically integrating Eq. 14 at 
24 local times, viz., t = 00, 01, 02, .... 23 hours. For an ionization chamber 


F(@, a) = 1. For counter telescope it can be evaluated by the formule 
given by Newell and Pressly.”° 


V. NUMERICAL EVALUATION OF THE NATURE OF THE SOLAR 
DAILY VARIATION OF MESON INTENSITY 


To compare the daily variation of meson intensity produced by the ani- 
sotropy of primary cosmic radiation with that observed experimentally, we 
must first examine some physical characteristics of the beam. The velocity 
of the ionized matter in the beam can be calculated from the time delay 
between the visible solar disturbances and the subsequent magnetic disturb- 
ances in the earth’s magnetic field. A delay of about 24 hours when there 
is a large magnetic storm following a solar flare, indicates a velocity of 
2,000 km./sec., while a delay from two to four days in the case of recurring 
magnetic disturbances suggest a velocity in the range of 1,000 km./sec. to 
500 km./sec. We assume therefore that the velocity of the ionized matter 
at different times varies between 2,000 km./sec. and 500 km./sec. 


The width of the beam is taken as ~6 x 10!2cm. With the synodic 
period of the sun as 27 days and the earth’s orbital radius as 1-5x 1013 cm., 
the earth will situate in the beam for about two days. The value is therefore 
consistent with inferences from geomagnetic data. 


At the present moment there is no method by which we can determine 
the strength of thefmagnetic field in the beam. Parker®! has recently shown 
that the value of the magnetic field in the beam is very insensitive to the 
density of the gas at the time of emission and with interplanetary density 
~ 10-*! g./cm.® it should be less than 1-4x10~ gauss. We assume that its 


value near the ecliptic varies between 2x 10-> gauss to 5x 10-* gauss. This 
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is necessary to produce an anisotropy compatible with experimental deter- 
mination of the solar daily variation at low latitudes which has an amplitude 
of 0-3 to 1-0%. 


With these values we show in Figs. 5 and 6 the expected daily variation 
of meson intensity at sea-level, i.e., atmospheric depth h = 1030 g./cm.* 
when the earth is just outside the beam and when the earth is within the beam 
at its centre respectively. For m(E, 1030) i(E), the values given by Naga- 
shima2? using the primary spectrum given by Neher as 


0-048 


i(E) = E23 (1 + 0-09 E3372 


are used. The amplitudes and the times of maxima of the Fourier compo- 
nents of the daily variation are given in Table I. It must be emphasized 
that these values depend upon i(E) dE, the differential energy spectrum of 
the primary radiation and m(E, 1030), the form of the multiplicity function. 
Recent evidence by Meyer and Simpson?? and Neher! suggests that the differ- 
ential energy spectrum varies during the course of years and a steeper spectrum 
than what is now assumed will give rise to larger amplitudes than those given 
in Table I. The precise value of the multiplicity function which is derived 
from the experimental determination of the latitude effect can be obtained for 
primary energies upto 15 Bev. For higher primary energies various authors 
have extrapolated it with different functions. Since at low latitudes primaries 
above 15 Bev. only are involved this extrapolation presents uncertainties. 
Hence the actual values at any time may be different from those given in 
Table I. However, they are not likely to differ much as to alter the picture 
presented therein. 


An important point that emerges from Figs. 5 and 6 is that while a posi- 
tive component of the trapped magnetic field within the beam gives rise to 
a diurnal variation at low latitudes with maximum shortly before noon, the 
reversal of this component will make the maximum in the daily variation 
shift to early morning hours. The daily variation in equatorial latitudes 
should therefore be expected to exhibit maximum either during night at about 
03 hours or during day at about 11 hours. Beams emitted at different periods 
may not have identical characteristics and the time of maximum may vary 
within an interval of about an hour. 


When the earth is at some distance from the beam, the walls of the beam 
will be inclined to the sun-earth line and the anisotropy in directions lying 
close to the walls of the beams is not accessible from the earth. The aniso- 
tropy of particles having energy near W.,,,, (0, 5) is therefore not observed 
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Fic. 5. The daily variation of meson intensity measured by an ionization chamber at sea- 
level in equatorial latitudes when the earth is just outside the beam of width ~ 6x 10!2cm. 





and the observed duration of the anisotropic region decreases. As the 
beam approaches the earth the accessible region increases, attains a maximum 
value within the beam and then decreases as the beam recedes from the earth. 
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Fic. 6. The daily variation of meson intensity measured by an ionization chamber at sea- 
level in equatorial latitudes when the earth is at the centre of the beam of width ~ 6 x 10'* cm, 


C The amplitudes of the variation which depends upon the effective duration 
1 of the anisotropy increases for a few days, reaches a maximum when the 








earth is at the centre of the beam and then decreases progressively. 
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The limit of the region of acceleration or deceleration given by Eqs. 5-12 
is expressed with reference to the walls of the beam. Their inclination to 
the sun-earth line introduces a correction to Ty ,. If we assume tentatively 
that the influence of the beam extends upto the time that it is within three 
days from the earth on either side of it, we have a maximum inclination of 
about 45° to the sun-earth line. The time of maximum would thus be dis- 
placed upto 3 hours on either side of the values indicated in Table I. Hence 


there should be a group of days when the time of maximum progressively 
shifts to later hours. 


A further factor to be considered is the inclination of the walls of the 
beam to the sun-earth line due to the radial emission of solar matter from the 
spinning sun. Chapman and Bartels*® have shown that the inclination would 
be by an angle of 12°, 24° and 48°, if the velocity of the matter is 2,000 km./sec., 
1,000 km./sec. and 500 km./sec. respectively. The time of maximum would 
thus be further advanced by 12/15, 24/15 and 48/15 hours respectively. 


It will now be realised that the daily variation observed on individual 
days will not always resemble the curves shown in Figs. 5 and 6, but these 
can still serve to represent the average characteristics of the daily variation. 
Analysis of Huancayo data by Sarabhai et a/.* showing that the broad features 
of the daily variation of meson intensity at Huancayo can be explained in 
terms of the two types of variations, one having maximum at 1100 or 1300 
hours and the other having a maximum at 0100 or 0300 hours can now be 
understood. The time of maximum on individual days will not always be 
sharply at the values given in Table I but will be grouped round these values. 
Observations by Sarabhai and Nerurkar* on the daily variation of meson 
intensity at Ahmedabad show similar characteristics. 


The 12 monthly mean daily variation of meson intensity exhibits a 22-year 
cycle of change as observed at Huancayo. [If this is due to change in the rela- 
tive frequency of the variations with day and night maxima from year to year, 
the average orientation of the trapped magnetic field within beams should 
also have a similar variation. The most prominent feature on the sun which 
exhibits a 22-year cycle is the magnetic field associated with active solar 
regions. Maunder?* has shown from a separate plot of sunspot areas for 
each hemisphere that there are comparatively large differences in the sun- 
spot numbers of the two hemispheres as well as their time of maxima. One 
hemisphere appears to be more active than the other and remains so in the 
next sunspot cycle even though the reversal of the polarity of the magnetic 
fields associated with active regions takes place. If there is similar difference 
in activity between the North and the South hemispheres of the sun for 
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corpuscular emission, and the coherent magnetic field trapped within the 
beam is in some way related to the local fields of active solar regions it would 


be possibleto have the 22-year cycle of change in the daily variation of 
meson intensity. 


The amplitude of the daily variation on individual days can be as high 
as 0-7% for meson intensity as shown in Table I. It may be higher than 0-7% 
if the energy spectrum of the primaries is steeper. The amplitude averaged 
over a group of days, however, depends upon the number of days on which 
the beams are effective as well as the types of beams. Taking this number 
to be less than half it would be possible to account for the average amplitude 
of 0-1 to 0-3% for ionization chambers. 


With a particular orientation of the magnetic field within the beam, the 
type of anisotropy and hence the daily variation is similar when the earth is 
on either side of it or within it. Each type of daily variation is thus expected 
to occur in groups of days. The recurrence tendency is a natural consequence 


of the persistence over several solar rotations of an active region responsible 
for a beam. 


From a careful examination of geomagnetic disturbances and_ their 
association with solar phenomena in 1952-53, a period of low sunspot acti- 
vity when the active regions were well defined and well observed, Pecker and 
Roberts?? have shown that there is a minimum of geomagnetic activity about 
three days after the CMP of an active region. There are also definite maxima 
symmetrically situated at about 3 days from the minimum. They suggest 
that the solar active regions possess a magnetic field which deflects the ionized 
matter emitted by the region so that there is a deficiency of emitted matter 
at solar longitudes radially along the active centre. The ionized matter flows 
along lines of the magnetic field. At some distance above the centre of acti- 
vity, the kinetic energy of the ions can be sufficient to overcome the magnetic 
field so that the guiding effect of the field of the active centre controls the 
ion flow only near the centre. The angle of aperture of this cone of avoidance 
may vary with the strength of the magnetic field of the active centre as well 
as that of the emission. With this picture of emission of ionized matter 
from the sun, we have in the plane of the ecliptic two beams separated on the 
average by about 5 to 6 days. 


If the component of the trapped magnetic field perpendicular to the 
plane of the ecliptic is the same in both the beams and the earth is situated 
between them, it is seen from Figs. 3 and 4 that the earth receives accelerated 
particles from one beam and decelerated particles from the other. The ampli- 
tude of the daily variation can then be about double the values given in 
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Table I. If however this component of the magnetic field is different in both 
the beams, the carth receives accelerated or decelerated particles from each 
of them. Both types of anisotropies are simultaneously present. The result- 
ant variation depends upon the distance between the two beams as well as 
the characteristics of them and in some cases it is possible to have a daily 
variation with two maxima. When the earth is on one side of both the beams, 
only the adjacent beam is effective and we have the resultant variation with 
only one maximum. It is thus possible to accommodate on the present 
model the existence of a daily variation with more than one maximum, such 
as is reported by Sarabhai and Nerurkar.! 


VI. CONCLUSION 


The creation of an anisotropy by the presence in the neighbourhood 
of the earth of beams of solar ionized matter with trapped magnetic fields 
derived from local fields related to active regions on the sun furnishes a 
plausible interpretation of several significant features of the daily variation 
of meson intensity. If one adopts the density and the velocity of ionized 
matter in the beams consistent with solar relationships of geomagnetic 
disturbances and assumes a strength of the trapped magnetic field com- 
patible with electrodynamical considerations it is possible to derive a daily 
variation of meson intensity which, at low latitudes, would have a maximum 
just before noon or in the early morning. The two types of daily variation 
would occur in groups of days which would have a tendency of 27-day 
recurrence. The daily variation could have an amplitude upto 1% on 
individual days. The amplitude as well as the time of maximum of each 
type of daily variation would be expected to alter during each group. There 
would also be a possibility of daily variation with two maxima on some 
occasions. 


The present interpretation does not deal with a permanent anisotropy 
of cosmic rays about which there is at present little experimental knowledge 
but which may nevertheless be present to a small degree. On the other hand 
it attempts to explain the observed features of a highly variable anisotropy in 
terms of a modulation of the primary cosmic ray intensity. Consequently 
the changes in the daily variation are linked with day-to-day changes in the 
daily mean intensity of cosmic radiation. A crucial test of the present inter- 
pretation would lie in the study of the correlated behaviour of the changes 
in the daily variation and in the mean intensity. This examination is in pro- 
gress, 
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SUMMARY 


The present status of knowledge concerning the daily variation of meson 
intensity is reviewed. A solar daily variation of a constant nature arising out 
of a permanent anisotropy in the cosmic radiation has not so far been estab- 
lished and even if it exists it may not exceed a few tenths of a per cent. How- 
ever a number of features of a solar daily variation of a highly variable 
character having on the average a time of maximum on individual days either 
in the early morning or near noon are known. An attempt is therefore made 
to interpret the occurrence of a variable anisotropy of cosmic rays to which 
are attributed the principal experimental observations. 


Beams of solar ionized matter with a frozen magnetic field derived from 
a solar dipole field have been earlier suggested to explain abnormal diurnal 
variation on magnetically disturbed days through an anisotropy in primary 
cosmic radiation created by them. The present paper considers that the 
variable anisotropy is due to such beams but the frozen magnetic field is 
derived from the magnetic field of the active solar regions and bears no pre- 
ferential orientation in respect to the solar dipole field. With different 
orientations of the magnetic field in the beam, the component perpendicular 
to the solar equatorial plane could be reversed in direction. Depending upon 
the direction of this component, the anisotropy in the primary radiation gives 
rise to a daily variation at equatorial stations with maximum either in early 
morning hours or near noon for a group of days. With the strength of this 
component of the magnetic field in the range 10-* to 10-* gauss and the velocity 
of the beam between 500 km./sec. and 2,000 km./sec., the daily variation of 
meson intensity has amplitudes from 0-3 to 1-5%. The simultaneous pre- 


sence of two such beams gives rise to double the amplitudes or a daily varia- 
tion with two maxima. 


Some further applications of this concept are discussed. 
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